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INTRODUCTION 

Latest observation of la type supernovas and CMB yields to following compo- 
sition of the Universe: 4% of baryonic matter, 20% of dark matter and 76% 
of dark energy. The term "dark matter" is implemented for the unknown 
matter, which has an ability for clustering, but was not yet detected in lab- 
conditions. The "dark energy" is the energy not been detected yet, but also 
unable for clustering, as common energy does. 

Most probably, the dark energy is responsible for the value of the cosmo- 
logical constant. Recent experiments indicate the value of the cosmological 
constant to be too small if originated only by vacuum energy of common 
matter. This brings some difficulties for developing of corresponding theo- 
ries. A few ways for interpretation of presents for the cosmological constant 
are known. Some of those are dynamical alternatives of dark energy. 

In this Thesis the attention is focused on f{R) theory. This kind of grav- 
itational theory is the result of the lagrangian generalization in the Hilbert- 
Einstein action. On the base of this theory a study will be undertaken on 
the problem of the effective cosmological constant, accelerated expansion of 
the Universe and the compactification (non-observation) of additional dimen- 
sions. 

Topicality of the subject. 

Multidimensionality of our Universe is one of the most intriguing assump- 
tion in modern physics. It follows naturally from theories unifying different 
fundamental interactions with gravity, e.g. M/string theory The idea 

has received a great deal of renewed attention over the last few years. How- 
ever, it also brings a row of additional questions. 

According to observations the internal space should be static or nearly 
static at least from the time of primordial nucleosynthesis, otherwise the 
fundamental physical constants would vary. This means that at the present 
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evolutionary stage of the Universe there are two possibilities: slow variation 
or compactification of internal space scale parameters [3] . 

in many recent studies the problem of extra dimensions stabilization was 
studied for so-called ADD (see e.g., Refs. P-[TT]). Under these approaches 
a massive scalar fields (gravitons or radions) of external space-time can be 
presented as conformal excitations. 

In above mentioned works it was assumed that multidimensional action to 
be linear with respect to curvature. Although as follows from string theory, 
the gravity action needs to be extended to nonlinear one. In order to investi- 
gate effects of nonlinearity, in this Thesis a multidimensional Lagrangian will 
be studied, having the form L = f{R), where f{R) is an arbitrary smooth 
function of the scalar curvature. 

Connection of the Thesis with scientific programs, 
plans and projects. 

The study under this Thesis was undertaken as part of the budget project 
of Odessa National Mechnikov University #519 "Compactification of space- 
time in quantum cosmology". 

The target and the task of the study. 

The task of this study is an investigation of diverse non-linear cosmological 
models with respect to condition of additional dimensions compactification 
and an analysis of additional dimensions effect on the evolution of the Uni- 
verse. 

For that, several surveys have to be implemented: an analysis of non- 
linearity influence on the effective potential in the equivalent theory; ca- 
pability of extrema for these potentials as result of additional dimensions 
non-observability condition. Stabilization and compactification of additional 
dimensions and their consistency with observations. 

Hence, the attention is accumulated on the possibility to achieve models, 
in which external space (our 4D space-time) behaves as being observed Uni- 
verse and internal space is stabilized and compactified on the Plank scales. 
These kind of models allow to build a bridge between observations, multi- 
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dimensionality and evolution of the Universe, proving by this the principal 
possibility for existence of additional dimensions. 

The scientific novelty of obtained results. 

The multidimensional models of type + for pure geometric action 
and of with forms are studied for stabilization and compactification of 
additional dimensions. 

It is shown, that for non-linearity of type in multidimensional case 
with forms, the stabilization of internal spaces is independent on signatures 
of internal space curvature, multidimensional and effective cosmological con- 
stants. Moreover, effective cosmological constant may satisfy to observed 
densities of dark energy by applying fine tuning. First time for this kind 
of model, the capability for stable compactification of additional dimensions 
and accelerated expansion of the Universe is shown. When in case of pure 
geometrical model of R~^ this two phenomena simultaneously are not obtain- 
able. 

The analysis of infiation is undertaken for linear, quadratic, quartic models 
with forms and for R^ + R^ model with pure gravitational action. 

A new type of inflation is found for + R^ model, called as bouncing 
inflation. For this type of inflation there is no need for original potential to 
have a minimum or to check the slowroll conditions. A necessary condition 
is the existence of the branching points. It is shown that this inflation takes 
place both in the Einstein and Brans-Dicke frames. 

Practical value of obtained results. 

The obtained results have fundamental and theoretical value, giving pos- 
sibility for assay of the Universe evolution, taking into account presence of 
additional dimensions for diverse non-linear models either pure gravitational 
or with forms. The parameters values obtained in the Thesis may become an 
object for further study and justification in experiments. 

Approbation of the Thesis results. 

The main results of the study, present in the Thesis, were reported on 
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CHAPTER 1 
LITERATURE OVERVIEW 

The currently observed accelerated expansion of the Universe suggests that 
cosmic flow dynamics is dominated by some unknown form of dark energy 
characterized by a large negative pressure [12]. Most recent observations of 
CMB and la supernovas indicates that Universe consists of: 76% of dark 
energy, 20% of dark matter and 4% of baryonic matter jlBHIS]- Many ap- 
proaches are known for distinguishing of dark matter and dark energy. Usu- 
ally it is assumed that dark energy despite of dark matter is incapable for 
clustering and satisfies for Strong Energy Condition |I7l[IB] • Due to dark en- 
ergy dominates over matter, it is responsible for late-time acceleration of the 
Universe (i.e. determines the value of effective cosmological constant), follow- 
ing after early-time acceleration as predicted with inflation paradigm ||T2H2I] • 
In between these to stages, the period of decelerated expansion is needed to 
provide eras of radiation domination for BBN and matter domination for 
formation of Universe structure 

As general approach, the standard Einstein's General Relativity (GR) is 
used as paradigm (foundation) for developing of new theories, called as Ex- 
tended Theories of Gravity (ETG). This allows to keep successful features of 
GR and apply extensions and corrections to Einstein's theory. Usually by 
adding higher order curvature invariants and/or minimally or non-minimally 
coupled scalar fields to the dynamics; these corrections emerge from the ef- 
fective action of quantum gravity [I21I22] • 

Among many, the simplest theory capable more or less to describe ex- 
perimental data is so-called ACDM model. It gives appropriate qualitative 
picture of the observed Universe, but do not explain the inflation. Also, in 
different models well known problem of the cosmological constant magnitude 
arises being higher in 120 orders of magnitude then observed one j23] . 
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As next step for modification of GR, one may apply Mach's principle, 
which states that the local inertial frame is determined by the average mo- 
tion of distant astronomical objects [24]. This means that the gravitational 
coupling could be determined by the distant distribution of matter, and it 
can be scale-dependent and related to some scalar field. As a consequence, 
the concept of "inertia" and the Equivalence Principle have to be revised. 
Brans-Dicke theory [25] constituted the first consistent and complete the- 
ory alternative to Einstein's GR. Brans-Dicke theory incorporates a variable 
gravitational coupling strength whose dynamics are governed by a scalar field 
non-minimally coupled to the geometry, which implements Mach's principle 
in the gravitational theory jT2l[25H27] . 

Nonlinear ETG models may arise either due to quantum fluctuations of 
matter fields including gravity [28] , or as a result of compactification of extra 
spatial dimensions [2Sj- Compared, e.g., to others higher-order gravity the- 
ories, f{R) theories are free of ghosts and of Ostrogradski instabilities [30]. 
Recently, it was realized that these models can also explain the late time ac- 
celeration of the Universe. This fact resulted in a new wave of papers devoted 
to this topic (see e.g., recent reviews II21EIIE2])- 

Another intriguing assumption in modern physics is the multidimensional- 
ity of our Universe. It follows from theories which unify different fundamental 
interactions with gravity, such as M or string theory |I],[2], and which have 
their most consistent formulation in spacetimes with more than four dimen- 
sions. Thus, multidimensional cosmological models have received a great deal 
of attention over the last years. 

Stabilization of additional dimensions near their present day values (dila- 
ton/geometrical moduli stabilization) is one of the main problems for any 
multidimensional theory because a dynamical behavior of the internal spaces 
results in a variation of the fundamental physical constants. Observations 
show that internal spaces should be static or nearly static at least from the 
time of recombination (in some papers arguments are given in favor of the 
assumption that variations of the fundamental constants are absent from 
the time of primordial nucleosynthesis [SB])- In other words, from this time 
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the compactification scale of the internal space should either be stabilized 
and trapped at the minimum of some effective potential, or it should be 
slowly varying (similar to the slowly varying cosmological constant in the 
quintessence scenario). In both cases, small fluctuations over stabilized or 
slowly varying compactification scales (conformal scales/geometrical moduli) 
are possible. 

Stabilization of extra dimensions (moduli stabilization) in models with 
large extra dimensions (ADD-type models) has been considered in a number 
of papers (see e.g., Refs. |lll[II]). In the corresponding approaches, a product 
topology of the (4 + D')— dimensional bulk spacetime was constructed from 
Einstein spaces with scale (warp) factors depending only on the coordinates 
of the external 4— dimensional component. As a consequence, the confor- 
mal excitations had the form of massive scalar fields living in the external 
spacetime. Within the framework of multidimensional cosmological models 
(MCM) such excitations were investigated in [SIEB] where they were called 
gravitational excitons. Later, since the ADD compactification approach these 
geometrical moduli excitations are known as radions [HE] . 

Most of the aforementioned papers are devoted to the stabilization of 
large extra dimensions in theories with a linear multidimensional gravita- 
tional action. String theory suggests that the usual linear Einstein-Hilbert 
action should be extended with higher order nonlinear curvature terms. In 
the papers [ j551[57] a simplified model is considered with multidimensional La- 
grangian of the form L = f{R)-, where f{R) is an arbitrary smooth function 
of the scalar curvature. Without connection to stabilization of the extra- 
dimensions, such models (4— dimensional as well as multidimensional ones) 
were considered e.g. in Refs. [2BH1B]- There, it was shown that the nonlinear 
models are equivalent to models with linear gravitational action plus a min- 
imally coupled scalar field with self-interaction potential. Similar approach 
was elaborated in Refs. [47] where the main attention was paid to a possi- 
bility of the late time acceleration of the Universe due to the nonlinearity of 
the model. 

The most simple, and, consequently, the most studied models are polyno- 
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mials of R: f{R) = Yllli^Q^nR^ [k > 1), e.g., quadratic R + R^ and quartic 
R-\-R^ ones. Active investigation of these models, which started in 80-th years 
of the last century [1BH5D], continues up to now fMlES]. Obviously, the cor- 
rection terms (to the Einstein action) with n > 1 give the main contribution 
in the case of large R, e.g., in the early stages of the Universe evolution. As 
it was shown first in [53] for the quadratic model, such modification of grav- 
ity results in early inflation. From the other hand, function f{R) may also 
contain negative degrees of R. For example, the simplest model is ^ + R^. 
In this case the correction term plays the main role for small R, e.g., at the 
late stage of the Universe evolution (see e.g. fSZllSl] and numerous references 
therein). Such modification of gravity may result in the late-time acceleration 
of our Universe [55]. Nonlinear models with polynomial as well as ^^^-type 
correction terms have also been generalized to the multidimensional case (see 

e.g., [szliaiisMi]). 

In this Thesis several models of /(^)-theory are analyzed for pure gravi- 
tational case in chapter 2, and with forms in chapter 3, as well as inflation 
possibility is studied in chapter 4. 
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CHAPTER 2 



ASYMPTOTICALLY AdS NONLINEAR 
GRAVITATIONAL MODELS 

In this chapter the non-hnear gravitational models with curvatures of type 
R"^ and + R^ are studied. It is shown that for particular values of pa- 
rameters, the stabilization and compactification of additional dimensions is 
achievable with negative constant curvature of internal space. In this case 
the 4— dimensional effective cosmological constant becomes negative. As re- 
sult homogenous and isotropic external time-space appears to be AdS4. The 
correlation between D— dimensional and 4— dimensional fundamental masses 
imposes restrictions on the parameters in the models. 



2.1. General setup 



Let us consider a D = (4+ D')— dimensional nonlinear pure gravitational 
theory with action functional 



S = Wjr I d^'x^WiR), (2.r 



2^1) J M 

where f{R) is an arbitrary smooth function with mass dimension 0{inn?) 
(m has the unit of mass) of a scalar curvature R = R[g] constructed from 
the D— dimensional metric gab (a, 6 = 1, . . . , /^). D' \s the number of extra 
dimensions and ti\) denotes the D— dimensional gravitational constant which 
is connected with the fundamental mass scale M^(4+£)/) and the surface area 
Sd-1 = 271^^^^^/"^ /T[{D — l)/2] of a unit sphere m. D — 1 dimensions by the 
relation E2 



4 = 2SD-i/Mll,%y (2.2) 

Before endowing the metric of the pure gravity theory ( 12.11 ) with explicit 
structure, let us recall that this ^—nonlinear theory is equivalent to a theory 
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which is hnear in another scalar curvature R but which contains an additional 
self-interacting scalar field. According to standard techniques [38-46. 50] , the 
corresponding i?— linear theory has the action functional: 

^ = ?V / d^'x^l [R[g] - ^«V,a0,6 - 2t/(0)] , (2.3) 

where 

f^n)^^:^e^^>0, (2.4) 
and where the self-interaction potential U (0) of the scalar field is given by 

= le-''^\R{(j))e^^-f(R{(j)))],B:= ^ ^ (2.5) 
2 ^ ^ ^ ^^^^^ ^{D-2){D-1) ^ ^ 

This scalar field (p carries the nonlinearity degrees of freedom in R of the 
original theory, and for brevity will be called the nonlinearity field. The 
metrics gab, gab of the two theories (12.11 ) and (12.31 ) are conformally connected 
by the relations 

gab = n%b=[f'{R)f^''~'^-9ab. (2.6) 
Next, let us assume that the D-dimensional bulk space-time M undergoes 
a spontaneous compactification to a warped product manifold 

M = Mq X Ml X . . . X Mn (2.7) 

with metric 

n 

g = gab{X)dX- dX' = ^(0) + ^ e'^'^^^g^'^ . (2.8) 

i=i 

The coordinates on the {Dq = do + 1)— dimensional manifold Mq (usually 
interpreted as our observable {Dq = 4)— dimensional Universe) are denoted 
by X and the corresponding metric by 

^(0) = g(^J{x)dx^'(g)dx'' . (2.9) 

For simplicity, the internal factor manifolds Mi are chosen as (i— dimensional 
Einstein spaces with metrics g^"^^ = gm,n,{yi)dy^' dy^\ so that the relations 

= >''9^w m„n, = l,...,6i, (2.10) 



R 



mitii 
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and 



R 



= AV, = R, (2.11) 

hold. The specific metric ansatz (12.81 ) leads to a scalar curvature R which 
depends only on the coordinates x of the external space: R[g] = R{x). Cor- 
respondingly, also the nonlinearity field depends on x only: = (f){x). 

Passing from the ^—nonlinear theory (12.11 ) to the equivalent i?— linear 
theory (12.31 ) the metric (12.81 ) undergoes the conformal transformation g ^ g 
[see relation (12.61 )] 

n 

g^n'-g^ (e-^*)'""-''fl := 5<°' + ^e^^-WflO (2.12) 



i=l 



with 



gV^ := ie^*f^^-'^ gi^ , /?-:=/?- + ^0. (2.13) 
2.2. Stabilization of internal dimensions 

The main subject of subsequent considerations will be the stabilization of 
the internal space components. A strong argument in favor of stabilized or 
almost stabilized internal space scale factors /3*(a:), at the present evolution 
stage of the Universe, is given by the intimate relation between variations 
of these scale factors and those of the fine-structure constant a The 



strong restrictions on a— variations in the currently observable part of the 
Universe jMHST] imply a correspondingly strong restriction on these scale 



factor variations [53]. For this reason, the derivation of criteria ensuring a 
freezing stabilization of the scale factors will be performed below. 

In Ref. [25] it was shown that for models with a warped product structure 
(12.81 ) of the bulk spacetime M and a minimally coupled scalar field living on 
this spacetime, the stabilization of the internal space components requires a 
simultaneous freezing of the scalar field. Here a similar situation with simul- 
taneous freezing stabilization of the scale factors P^{x) and the nonlinearity 
field (l){x) is expected. According to (12.131 ). this will also imply a stabihzation 
of the scale factors ^*(x) of the original nonlinear model. 
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In general, the model will allow for several stable scale factor configurations 
(minima in the landscape over the space of volume moduli). Let us choose 
one of them§, denote the corresponding scale factors as /^g, and work further 
on with the deviations 

^\x) = ^\x) - (31, (2.14) 

as the dynamical fields. After dimensional reduction of the action functional 
(12.31 ) one passes from the intermediate Brans-Dicke frame to the Einstein 
frame via a conformal transformation 

/ n \ -2/P0-2) 
<i = «^9<2= n«*'' (2.15) 



with respect to the scale factor deviations ^^{x) [SSIIS^- As result the follow- 
ing action is achieved 



^(°)'^^5^(/.5,(^-2[/e//|, (2.16) 



Mo 



which contains the scale factor offsets /^g through the total internal space 
volume 

n „ n 

Vd'^ViXVo = \{ dS\l b«| X n e^'^o (2.17) 

^=l ^^^^ i=\ 

in the definition of the effective gravitational constant of the dimension- 
ally reduced theory 

4-7r 

«^(z?o^4) = = Stt/MI =^ M! = —Vn'M'^+%,y (2.18) 

Obviously, at the present evolution stage of the Universe, the internal space 
components should have a total volume which would yield a four-dimensional 
mass scale of order of the Planck mass M(4) = Mpi. The tensor components 
of the midisuperspace metric (target space metric on r^) reads: Gij = diSij + 



^Although the toy model ansatz (|2.1I) is highly oversimplified and far from a realistic model, one can 
roughly think of the chosen minimum, e.g., as that one which is expected to correspond to current evolution 
stage of our observable Universe. 



18 



didj/ {Dq — 2), where i, j = {1, . . . , n), see pSlE2]- The effective potential has 
the explicit form 



Ueff{f3,<P) = 

where abbreviated 



n 



-Dn-2 



-2/3* 



+ ^(0) 



i=l 



(2.19) 



;2.20) 



R, := i?,exp(-2/3a. 

A freezing stabilization of the internal spaces will be achieved if the ef- 
fective potential has at least one minimum with respect to the fields (5'^{x). 
Assuming, without loss of generality, that one of the minima is located at 

= /?Q ^ /3* = 0, the extremum condition reads: 



dU, 



eff 



= 



Ri = 



d, 



/3=0 



Dr 



Y.Rj + 2U{ct>) 



:2.2i) 



From its structure (a constant on the l.h.s. and a dynamical function of 
on the r.h.s) it follows that a stabilization of the internal space scale factors 
can only occur when the nonlinearity field 0(x) is stabilized as well. In the 
freezing scenario this will require a minimum with respect to 0: 



dU{(j)) 



^ 



dU, 



eff 







:2.22: 



Hence, a stabilization problem is arrived, some of whose general aspects have 
been analyzed already in Refs. ||31HSZ]- For brevity let us summarize the 
corresponding essentials as they will be needed for more detailed discussions 
in the next consideration. 



1. Eq. ( 12.211 ) implies that the scalar curvatures Ri and with them the com- 
pactification scales e^o [gee relation ( 12.201 )] of the internal space compo- 
nents are finely tuned 



Ri 
d, 



Rj 

d-i 



i,j = l,...,n 



[2.23) 
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2. The masses of the normal mode excitations of the internal space scale 
factors (gravitational excitons/radions) and of the nonlinearity field 
near the minimum position are given as [33]: 



mi 



ml = -^^U{<J^o) = -2^ > , (2.24) 



2 d'Ui^p) 
ml := - 



> 



(2.25) 



3. The value of the effective potential at the minimum plays the role of an 
effective 4D cosmological constant of the external (our) spacetime Mq: 



Ae// := U, 



eff 



;2.26) 



4. Relation (12.261 ) implies 



signAe// = signt/(0o) = signi?, 



;2.27) 



Together with condition (12.241 ) this shows that in a pure geometrical 
model stable configurations can only exist for internal spaces with neg- 
ative curvature^: 



R, < 



.,n). (2.28) 

Additionally, the effective cosmological constant Ag// as well as the min- 
imum of the potential t/(0) should be negative too: 



Aeff < 0, 



U{cj)o) < 



;2.29) 



Plugging the potential U{(f)) from Eq. (12.51 ) into the minimum conditions 



^Negative constant curvature spaces Mi are compact if they have a quotient structure: Mi = H'^' /Ti, 
where H'^^ and are hyperbohc spaces and their discrete isometry group, respectively. 
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(E22D, (E25D yields with the help of d^R = Af'/f" the conditions 



dU 



A 



2(D-2] 



if 



i.-D/{D-2) 



h 



= 0, 



h := Df-2Rf, 



M0o) = o, 



;2.30) 



= W^-''^^^[d^R+{A-B)R\^^ 



1 



(/ 



>o, 



where the last inequality can be reshaped into the suitable form 

f"dRh\^„ = f" [{D-2)f'-2Rf"]^^>0. 
Furthermore, from Eq. (12.301 ) follows 



(2.31) 



(2.32) 



^(0o) = ^(/O"^i^(0o 



so that (12.291 ) leads to the additional restriction 

i^(0o) < 



2.33) 



(2.34) 



at the extremum. 

Thus, to avoid the effective four-dimensional fundamental constant vari- 
ation, it is necessary to provide the mechanism of the internal spaces stabi- 
lization. In these models, the scale factors of the internal spaces play the role 
of additional scalar fields (geometrical moduli/gravexcitons [34]). To achieve 
their stabilization, an effective potential should have minima with respect to 
all scalar fields (gravexcitons and scalaron). Previous analysis (see e.g. [35]) 
shows that for a model of the form (12.31 ) the stabilization is possible only for 
the case of negative minimum of the potential U (0). According to the defini- 
tion (12.41 ). the positive branch of f'{R) > will be considered. Although the 
negative f'<0 branch can be considered as well (see e.g. Refs. P7l[5Ui[55] ). 
However, negative values of f'{R) result in negative effective gravitational 
"constant" Gg// = i^jj/ f- Thus /' should be positive for the graviton to 
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carry positive kinetic energy (see e.g., |I2])- From action ( 12.31 ) the equation 
of motion of scalaron field can be obtained in the form: 

dU 

□c/>-^ = 0. (2.35) 

0(p 

Now, let us analyze the internal space stabilization conditions (12.231 ) - 
(12.291 ) and (12.301 ) - (12.341 ) on their compatibility with particular scalar curva- 
ture nonlinearity f{R)- 



2.3. R-^ model 



Recently it has been shown in Refs. [7DH7B] that cosmological models with 
a nonlinear scalar curvature term of the type can provide a possible 
explanation of the observed late-time acceleration of our Universe within a 
pure gravity setup. The equivalent linearized model contains an effective 
potential with a positive branch which can simulate a transient inflation- 
like behavior in the sense of an effective dark energy. The corresponding 
considerations have been performed mainly in four dimensions. Here these 
analyses are extended to higher dimensional models, assuming that the scalar 
curvature nonlinearity is of the same form in all dimensions. Starting from a 
nonlinear coupling of the type: 

f{R) = R- li/R, /i>0, (2.36) 

in front of the term, the minus sign is chosen, because otherwise the 

potential U (0) will have no extremum. 

With the help of definition (12.41 ). the scalar curvature R can be expressed in 
terms of the nonlinearity field (p and obtain two real- valued solution branches 

R± = (e^'^ - 1) "^^^ , =^ > (2.37) 

of the quadratic equation f'{R) = e"^^. The corresponding potentials 

U±{(f)) = ±^e-^^ \/e^<^ - 1 (2.38) 

®A discussion of pro and contra of a higher dimensional origin of R^^ terms can be found in Ref. [29] . 
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have extrema for curvatures [see Eq. ( 12.301 )] 



'D + 2 
D-2 



e^-^o = ^ = > 1 for D>3 (2.39) 

and take for these curvatures the values 

The stabihty defining second derivatives [Eq. (12.311 )] at the extrema (12.391 ). 



show that only the negative curvature branch R yields a minimum with 
stable internal space components. The positive branch has a maximum with 
U+{(f)o) > 0. According to (12.271 ) it can provide an effective dark energy 
contribution with Ag// > 0, but due to its tachyonic behavior with d'^U (0o) < 
it cannot give stably frozen internal dimensions. This means that the 
simplest extension of the four-dimensional purely geometrical setup of 
Refs. fTDHTH] to higher dimensions is incompatible with a freezing stabilization 
of the extra dimensions. A possible circumvention of this behavior could 
consist in the existence of different nonlinearity types fi{Ri) in different factor 
spaces Mi so that their dynamics can decouple one from the other. This could 
allow for a freezing of the scale factors of the internal spaces even in the case of 
a late-time acceleration with Ag// > 0. Another circumvention could consist 
in a mechanism which prevents the dynamics of the internal spaces from 
causing strong variations of the fine-structure constant a. The question of 
whether one of these schemes could work within a physically realistic setup 
remains to be clarified. 

Finally, in the minimum of the effective potential t/e//(0? Z^*)' which is 
provided by the negative curvature branch R{(f)), one finds excitation masses 
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for the gravexcitons/radions and the nonhnearity field (see Eqs. ( 12.241 ). ( 12.401 ) 
and ( ICTD ) of order 

mi = . . . = rrin ~ ~ /i^^^ . (2.42) 



For the four-dimensional effective cosmological constant Ag// defined in ( 12.261 ) 
one obtains in accordance with Eq. ( 12.401 ) Ae// ~ —-s/J^- 

Summary 

As stability condition the existence of a minimum of the effective po- 
tential of the dimensionally reduced theory is assumed, so that a late-time 
attractor of the system could be expected with freezing stabilization of the 
extra-dimensional scale factors and the nonlinearity field. It was shown in 
Refs. [EZIISB], that for purely geometrical setups this is only possible for neg- 
ative scalar curvatures, (^ < 0), independently of the concrete form of the 
function f{R). 

Four-dimensional purely gravitational models with curvature contri- 
butions have been proposed recently as possible explanation of the observed 
late-time acceleration (dark energy) of the Universe [7DH7B] • It is shown that 
higher dimensional models with the same scalar curvature nonlinearity 
reproduce (after dimensional reduction) the two solution branches of the four- 
dimensional models. But due to their oversimplified structure these models 
cannot simultaneously provide a late-time acceleration of the external four- 
dimensional spacetime and a stabilization of the internal space. A late-time 
acceleration is only possible for one of the solution branches — for that which 
yields a positive maximum of the potential U (0) of the nonlinearity field. A 
stabilization of the internal spaces requires a negative minimum of U{(J)) as 
it can be induced by the other solution branch. 

2.4. R^ + model 

In this section let us analyze a model with curvature-quadratic and curvature- 
quartic correction terms of the type 

f{R) = R + aR^ + -fR^ - 2Ad . (2.43) 
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Figure 2.1: The surfaces /' = as a functions R = R{a,'^) for the model fl2.43p . 



According to eq. ( 12.41 ) 



1 + 2aR + 47i?^ 



;2.44) 



The definition ( 12.41 ) /' = eyL^{A(j)) clearly indicates that the positive branch 
/' > is chosen. For the given model ( 12.431 ). the surfaces /' = as a functions 
R = R{a, 7) are given in Fig. 12.11 As it easily follows from Eq. ( 12.441 ). points 
where all three values R^ 7 and a are positive correspond to the region /' > 0. 
Thus, this picture shows that there is one simply connected region /' > 
and two disconnected regions /' < 0. 

Eq. ( 12.441 ) can be rewritten equivalently in the form: 



a 



1 



:2.45) 



-00 < < +00 
-l< X < +00 . 



;2.46) 



Eq. ( 12.451 ) has three solutions ^^1,2,35 where one or three of them are real- 
valued. Let 
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The sign of the discriminant 

Q ■- r^ + q^ 
defines the number of real solutions: 



;2.48) 



g>o 

Q = 

Q<0 



^Rl = 0, ^i?2,3 ^ , 

^R^ = Vi, i?i = i?2 , 
^i?, = . 



(2.49) 



Physical scalar curvatures correspond to real solutions Ri{X). It is the most 
convenient to consider Ri = Ri{X) as solution family depending on the two 
additional parameters (a, 7): sign{a) = sign{'y) =4> Q > 0, sign{a) ^ 
sign{'j) =^ Q ^0. 

For Q > the single real solution ^1 is given as 



r + Q 



1/2 



1/3 



+ 



r-Q 



1/2 



nl/3 



Zl + ^2 



:2.50) 



where zi^2 is defined in the form: 

cosh(6') = 



2 2 
p = r 



Q = -q\ 



r 



Taking into account eq. ( 12.471 ). the function X reads 



X(^) = 87V^cosh(^). 
The three real solutions ^1^2,3 (-^) for Q < are given as 



(2.51) 



:2.52) 



Rl = Si + S2, 

R2 = ^(-l + ^V3)si + ^(-l-zV3) 



3 Si + e * 3 



i^3 = ^(-l-^V3)si + ^(-l + zV3) 



S2 



S2 



= e ' 3 + 3 S2, 



:2.53) 



where one can fix the Riemann sheet of Q^^^ by setting in the definitions of 

Sl,2 



Si 2 



r±i\Q\ 



1/2 



1/3 



(2.54) 
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A simple Mathematica calculation gives for Vieta's relations from (12.531 ) 



R1R2R3 = sl + sl = 2r. (2.55) 



In order to work with explicitly real-valued Ri let us rewrite si^2 from (12.541 ) 
as follows 



si,2 = K^'e 



\b\' = r'+\Q\=r'-Q = -q' 



and get via (12.531 ) 



R^ = si + S2 = 2|6|i/3cos(i?/3), (2.57) 
= e*fsi + e"*^S2 = 2|6|^/3 COS (7^/3 + 27r/3), 
= e-*^si + e*^S2 = 2|6|^/3 cos (7^/3 -27r/3) 



or 



Rk = 2|6|i/3cos 



i9 + 27r/c^ 



= 2v^cos (^^^^^ , /c = -1, 0, 1. (2.58) 

In order to understand the qualitative behavior of these three real-valued 
solutions as part of the global solution picture let us first note that, according 
to (12.451 ). one may interpret X as single- valued function 



X{R) = 4-fR^ + 2aR (2.59) 

and look what is happening when (a, 7) are changed. Obviously, the inverse 
function R{X) has three real- valued branches when X{R) is not a monotonic 
function but instead has a minimum and a maximum, i.e. when 

dnX ■- X' = 127i?2 + 2a = =^ R^ = -— (2.60) 

67 
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has two real solutions R± = zty^— a/ (67) and corresponding extremac 

X{R±) = ^aR± . (2.61) 

It should hold sign (a) 7^ sign (7) in this case, so that one find 

7 > 0, a < : Xmax = X{R^), X^^n = X{R+) 
7 < 0, a > : X^ax = X{R+), X^mn = X{R^) . 

(2.62) 

The transition from the three-real-solution regime to the one-real-solution 
regime occurs when maximum and minimum coalesce at the infiection point 

i?+ = i?_ = 0=^a = 0, 77^0. (2.63) 

(The non-degenerate case 7 7^ is considered here. Models with 7 = 
are degenerated ones and are characterized by quadratic scalar curvature 
terms only.) Due to — 1 < X < +00 the limit X — > +00 where in leading 
approximation may be considered: 

A-iR^ « X ^ +00 (2.64) 

so that 

R{X — > 00) — > sign (7) X 00 . (2.65) 

Leaving the restriction X > —1 for a moment aside, it was found that for 
product q;7 < there exist three real solution branches 7?.i^2,3: 

7 > : -00 <ni< R , -oo<X<Xr 

R- < ^2 < Xjnax ^ X > Xj 

R+ <'J^3< +00, Xjnin < X < +00, 

7 < : -00 <lZi< R , +00 > X > Xmin, 

R- < ^2 < R+, Xmin < X < Xjnaxi 
R+<'Tli< +00, Xmax > X > -00. 

(2.66) 



^m,axi 
^mini 



Hi is worth of noting that f"{R±) = X'{R±) = 0. 
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It remains for each of these branches to check which of the solutions Rk 
from (12.581 ) can be fitted into this scheme. Finally, one will have to set the 
additional restriction X > — 1 on the whole picture. 

To define the conditions for minima of the effective potential U^ff, first the 
extremum positions of the potential is to be obtained. The extremum 
condition ( 12.301 ) for the given particular model ( 12.431 ) reads: 



2.5. Case D = 8 - analytical solution 

In this part let us investigate the case of positive Q(0) that is equivalent 
to the condition 

Q{(j)) > =^ signa = sign7. (2.68) 

The case sign a ^ sign 7 that corresponds to different signatures of the 
discriminant Q will be considered in chapter 4. 

To define the conditions for minima of the effective potential 11^//, first 
let us obtain the extremum positions of the potential U{(J)). The extremum 
condition ( 12.301 ) for the given particular model ( 12.431 ) reads: 

^(0)i7 (y - 4) + ^(0)i« (y - 2) + ^(0)1 (f - 1) - = ' (2.69) 

where subscript " 1" indicates seeking the extremum positions for the solution 
( 12.501 ). Eq. ( 12.691 ) clearly shows that D = 8 is the critical dimension for the 
model ( 12.431 ) in full agreement with the result of the appendix 1 (see ( IA1.4P ). 



In what follows, that this critical case is investigated here. For D = 8 eq. 
( 12.691 ) is reduced to a quadratic one 



^(0)1 + - — = ; Ag = A^=8 (2.70) 

^CY CY 

with the following two roots: 




^lo)'i = -^±^/l-l . (2.71- 
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These roots are real if parameters a and Ag satisfy the following condition: 

If sign (a) = sign (Ag), then condition (12.721 ) is automatically executed, else 

l^sl < ^ , sign («) 7^ sign (Ag) . (2.73) 

To insure that roots (12.711 ) correspond to a minimum value of U {(/)), they 
should satisfy the condition (12.321 ): 

f" [{D - 2)/ - 2Rf"]^^ > ^ f"[3 + 4aR]^, > , (2.74) 

where 

/' = 2a + 12-fR^ . (2.75) 

Because for Q > eq. (12.501 ) is the single real solution of the cubic eq. 
(12.691 ). then R = Ri{(f)) is a monotonic function of Thus, the derivative 
dfj)Ri = Af'/f" does not change its sign. Keeping in mind that the f > 
branch is considered, the function Ri{(J)) is a monotone increasing one for 
f" > 0. As apparent form eq. (12.501 ). for increasing one should take 
7 > 0. In a similar manner, the function Ri{4>) is a monotone decreasing 
one for /",7 < 0. Thus, for the minimum position ^(o)i, inequality (12.741 ) 
leads to the following conditions (reminding that according to eq. (12.341 ) 
the minimum position ^(o)i should be negative and according to eq. (12.681 ) 
sign a = sign 7): 

I. r,7,«>0 : 

3 + 4«^fo)i > ^ |i?fo) J < 1^ . (2.76) 

II. r , 7, « < : 

3 + 4aRf,^, < ^ -|i?fo)J > ^ . (2.77) 



Obviously, inequality (12.771 ) is impossible and one arrives to the conclusion 
that the minimum of the effective potential f/e// is absent if sign a = sign 7 = 
-1. 
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Additionally, it can be easily seen that in the case 



sign a = sign 7 = signA^^i = +1 



;2.78) 



the effective potential Ue/f has no minima also. This statement follows from 
the form of the potential U{(f)) for the model (12.431 ). According to eq. (12.51 ). 
[/(0) reads: 

U{(j)) = (l/2)e-^'^ {aR^ + 87/?^ + 2Ai^) . (2.79) 

Thus, this potential is always positive for parameters satisfying (12.781 ) and 
one arrives to the contradiction with the minimum condition (12.291 ). There- 
fore, the investigation carried above indicates that the internal space stable 
compactification is possible only if the parameters satisfy the following sign 
relation: 

a > 0,7 > 0, Ag < 0. (2.80) 



Let us investigate this case in more detail. For this choice of signs of the 
parameters, it can be easily seen that both extremum values R^^^i from eq. 
satisfy the condition (ESI): R^^^ < 0. However, the expression 



>(±) ^ _ 

'(0)1 



l + 2aR^^^^+4-fR 



,(±)3 _ 

'(0)1 - 



1 9 

-±\ - - 16a|A8|-47 



R 



(±) 

(0)1 



:2.8i) 



shows that only R^^q^i can belong to /' > branch. To make /' 
positive, parameter 7 should satisfy the condition 



+ 



7 < 



16a|Ag 



R 



(+) 

(0)1 



(2.82) 



As apparent from this equation, parameter 7 remains positive if Ag belongs 
to the interval 



A« G 



1 

Sa 



,0 



;2.83) 



For this values of Ag, the condition (12.731 ) is automatically satisfied. It should 



be noted, that for positive a and negative R\^q^i the condition (12.761 ) is also 
satisfied. Taking into account the interval (12.831 ). the corresponding allowed 
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interval for 7 readc 



7G 



\ 



0, 



\ 



R 



(+) 

(0)1 



(2.84) 



Thus, for any positive value of a, Eqs. (12.831 ) and (12.841 ) define allowed 
intervals for parameters Ag and 7 which ensure the existence of a global 
minimum of the effective potential U^ff. Here, the required stable compact- 
ification of the internal space is obtained. The position of the minimum 
(/3q, 0o) and its value can be easily found (via the root R^^q^^) with the help of 
Eq. (12.441 ). The Figs.(I22D, (ESj) demonstrate such minimum for a particular 
choice of the parameters: a = 1,7 = 1,A8 = —0.1. Moreover, the limit 
Ag — > corresponds to r!(^^i — > which results in the decompactification of 
the internal space /^q — > 00. 



Summary 

In this section the model with curvature-quadratic and curvature-quartic 
correction terms of the type (12.431 ) was analyzed and it was shown that the 
stable compactification of the internal space takes place for the sign relation 
(12.801 ). Moreover, the parameters of the model should belong to the allowed 
intervals (regions of stability) (12.831 ) and (12.841 ). The former one can be 
rewritten in the form 

Ag = ^, ^£(-1,0). (2.85) 
8a 

Thus, for the root R^^^^i and parameter 7 one obtains respectively 



^(0)1 - 



r/ = i (-3 + ^9 + 8^) < (2.86) 



and g 

7 = ^, Ce(0,l). (2.87) 

Eq. ( ITM ) shows that r[q^^ G (-^.O)- 

^Similar interval for the allowed values of 7 was also found in [31, for the curvature-quartic model. 
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It is of interest to estimate the masses of the gravitational excitons (12.241 ) 



and of the nonlinearity field (f) (12.251 ) as well as the effective cosmological con- 



stant (I22SD. From Eqs. (^M)-(^M) follows that r\^.\ - Ag - -a'^ , 7 



10)1 

=^ /'(0o) ~ ^(1) ) /"(0o) ~ 5 U{(J)q) ~ —a^^. Then, the corresponding 
estimates read: 

- Aeff ^ ml ^ ml ^ . (2.88) 
From other hand (see Eqs. ([220]) and (E^l)) 



C/(0o) - exp(-2/3i) = 6(o2^ . (2.89) 

So, if the scale factor of the stabilized internal space is of the order of the 
Fermi length: 6(0)1 ^ Lp ^ lO^^^cm, then a ~ L^p and for the effective 
cosmological constant and masses one obtains: —^eff ~ ^1 ~ itt^I ~ ITeV^. 

The internal space stable compactification analysis was performed in the 
case Q{4>) > ^ sign a = sign 7. In forthcoming chapters the extend of this 
investigation to the case of negative Q{(j)) is performed, where the function 
R{(j)) has three real-valued branches. 

For the visualization reasons, let us consider the model with one internal 
space and critical dimension D = 8 (as usual, for the external spacetime 
Dq = 4). Then, the effective potential (12.191 ) reads: 



[2.90) 



To draw this effective potential, let us define ^1 via ?7(0o) in Eq. (12.241 ). In 
its turn, U{(j)o) is defined in Eq. ([23SD where R{(j)o) = r\qI and / (r^^^ 



(0)1 "^^^ J V (0)1 

can be found from (12.811 ). In Figs. ([2^, (12.31 ) the generic form of the Ue/f 
is illustrated by a model with parameters a = 1,7 = 1, Ag = —0.1 from the 
stability regions (ES) and ( ITM ). 
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Figure 2.2: Typical contour plot of the effective potential Ueff{f3^,4>) given in Eq. f l2.90p 
with parameters a = 1, 7 = 1, Ag = —0.1. f/e// reaches the global minimum at = 0, ~ 
-2.45). 




Figure 2.3: Typical form of the effective potential Ueff{f3^,(p) given in Eq. fl2.90p with 
parameters a = 1,7 = 1, Ag = —0.1. 
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CHAPTER 3 

MULTIDIMENSIONAL NONLINEAR MODELS 

WITH FORMS 

In the present chapter, the nonhnear gravitational multidimensional cosmo- 
logical model is considered with action of the type R+l/R with form-fields as 
a matter source, including a bare cosmological term as an additional param- 
eter of the theory. It is assumed that the corresponding higher-dimensional 
space-time manifold undergos a spontaneous compactification to a manifold 
with warped product structure of the external and internal spaces. Each of 
spaces has its own scale factor. A model without form-fields and bare cos- 
mological constant was considered in previous chapter, where the internal 
space freezing stabilization was achieved due to negative minimum of the 
effective potential. Thus, such model is asymptotically AdS without accel- 
erating behavior of our Universe. It is well known that inclusion of usual 
matter can uplift potential to the positive values [58]. One of the main task 
of present investigations is to observe such uplifting due to the form-fields. 
Indeed, it is demonstrated that for certain parameter regions the late-time 
acceleration scenario in present model becomes reachable. However, it is not 
simple uplifting of the negative minimum of the theory ( 12.11 ) to the positive 
values. The presence of the form-fields results in much more rich structure 
of the effective potential then for ( 12.11 ). Here, additional branches are ob- 
tained with extremum points, and one of such extremum corresponds to the 
positive minimum of the effective potential. This minimum plays the role 
of the positive cosmological constant. With the corresponding fine tuning 
of the parameters, it can provide the late-time accelerating expansion of the 
Universe. Moreover, it is shown that for this branch of the effective potential 
there is also a saddle point. Thus, domain walls are obtained, which separate 
regions with different vacua. It is demonstrated that these domain walls do 
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not undergo inflation because the effective potential is not flat enough around 
the saddle point. 

It is also worth of noting that the effective potential in this reduced model 
has a branchpoint. It gives very interesting possibility to investigate tran- 
sitions from one branch to another by analogy with catastrophe theory or 
similar to phase transitions in statistical theory. This idea needs more detail 
investigation. 

The chapter is structured as follows. First, a brief description of multi- 
dimensional models with scalar curvature nonlinearity f{R) and the form- 
flelds as a matter source is given. Then, dimensional reduction is per- 
formed and effective four-dimensional action with effective potential is ob- 
tained. General formulas from this section are applied to the speciflc model 
f{R) = R — ji/ R — Then minimum conditions of the effective potential 
are obtained. These conditions are analyzed for the cases of zero and posi- 
tive effective cosmological constants respectively. Furthermore, it is demon- 
strated that the positive minimum of the effective potential plays the role of 
the positive cosmological constant and can provide the late-time accelerating 
expansion. Additionally, this minimum is accompanied by a saddle point. It 
results in non-inflating domain walls in the Universe. 

3.1. General equations 

Let us consider a D = [Dq + D')— dimensional nonlinear gravitational 
theory with action functional 



where f{R) is an arbitrary smooth function of a scalar curvature R := R[g\ 
constructed from the D— dimensional metric gab (a, 6 = 1, . . . , i?). D' \s the 
number of extra dimensions. t^\) denotes the Z)— dimensional gravitational 



S 




(3.1) 
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constant. In action (13.11 ). a form field (flux) F has block-orthogonal struc- 
ture consisting of n blocks. Each of these blocks is described by its own an- 
tisymmetric tensor field F^'^\i = 1, . . . , n) of rank di (dj-form field strength). 
Additionally, let us assume that for the sum of the ranks holds Yl^^i di = D'. 

Following Refs. ||37l[79j, it can be shown that the nonlinear gravitational 
theory (13.11 ) is equivalent to a linear theory R = R[g] with conformally trans- 
formed metric 

gab = n%,= [f'{R)f^''-'^-g,, (3.2) 

and an additional minimal scalar field = ln[f'{R)]/A coupled with fluxes. 
The scalar field is the result and the carrier of the curvature nonlinearity of 
the original theory. Thus, for brevity, referring to the field as nonlinearity 
scalar field, a self-interaction potential U (0) of the scalar field reads 



/7(0) = ie-^^ [ i?(0)e^^ - / (i?(0))] , (3.3) 



where 



D-2 D , , 

A=x- , B= ^ 3.4) 

Furthermore, let us assume that the multidimensional space-time manifold 
undergoes a spontaneous compactification 

M — > M = Mq X Ml X . . . X Mn (3.5) 

in accordance with the block-orthogonal structure of the field strength F, 
and that the form fields F^'^\ each nested in its own dimensional factor 
space Mj (z = 1, . . . , n), respect a generalized Freund-Rubin ansatz [80]. Here, 
{Dq = 4)-dimensional space-time Mq is treated as our external Universe with 
metric g^^\x). 

This allows performing of a dimensional reduction of the model along the 
hues of Refs. j3S-13SlE2]. The factor spaces Mi are then Einstein spaces 
with metrics 5^^*^ = e^^^'^^^'j'^^^ which depend only through the warp factors 
ai{x) := e^*^^-* on the coordinates x of the external space-time Mq. For the 
corresponding scalar curvatures holds R [7^^*^] = X^di = Vi (in the case of the 
constant curvature spaces A* = ki{di — 1), ki = 0, ±1). The warped product 
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of Einstein spaces leads to a scalar curvature R which depends only on the 
coordinate x of the Do~ dimensional external space-time Mq: R[g] = R{x). 
This implies that the nonlinearity field is also a function only of x: (f) = (f){x). 
Additionally, it can be easily seen [27] that the generalized Freund-Rubin 
ansatz results in the following expression for the form-fields: (F*^*))^ = fi /^h^' 
where = const. 

In general, the model will allow for several stable scale factor configurations 
(minima in the landscape over the space of volume moduli). Let us choose 
one of them (which is expected to correspond to current evolution stage of 
our observable Universe), denote the corresponding scale factors as /^g, and 
work further on with the deviations ^^{x) = (5\x) — (3q. 

Without loss of generality§, one shall consider a model with only one 
(ii-dimensional internal space. After dimensional reduction and subsequent 
conformal transformation to the Einstein frame the action functional (13.11 ) 
readq^ 

S = d''^x^l\^\[R[g^'^^] --g^'^^^^d.ipd^if 

- ^(O)M.5,,05.0-2[/,;;(v9,0)} , (3.6) 

where (p := —y/di{D — 2)/[Dq — 2) and := i^\,/Vd^ denotes the Dq— 
dimensional (four-dimensional) gravitational constant. V^^ ~ exp((ii/3Q) is 
the volume of the internal space at the present time. 

A stable compactification of the internal space Mi is ensured when its 
scale factor (f is frozen at the minimum of the effective potential 



Uefj = e'^ 



■ii?ie"^ + [/(0) + /ie^V'^i^ 



(3.7) 



where R\ := riexp(— 2/3q) defines the curvature of the internal space at 
the present time and contribution of the form-field into the effective action 

^Thc difference between a general model with n > 1 internal spaces and the particular one with n — 1 
consists in an additional diagonalization of the geometrical moduli excitations. 

^The equivalency between original higher dimensional and effective dimensionally reduced models was 
investigated in a number of papers (see e.g. [HIHHS])- The origin of this equivalence results from high 
symmetry of considered models (i.e. because of specific metric ansatz which is defined on the manifold 
consisting of direct product of the Einstein spaces). 
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is described by h := nj^ exp {—2di(3Q) > 0. For brevity let us introduce 
notations 



a : = 2 \ I ——^ -T , b : = 2. ^ 



di{D-2)' • y (D-2)(Do-2) 

2di- D , , 

c : = = . 3.8 

V(D-l)(D-2) 



3.2. R model with forms 



In this section the conditions of the compactification are analyzed for a 
model with 

f{R) = R-^-2An. (3.9) 
Then from the relation f'{R) = exp {Acj)) one obtains 



R = g(^^A^ _ ^-^ ^ g = ±1 , s = sign(^) . (3.10) 

Thus, the ranges of variation of (j) are (j) G (— oo,0) for /i < (s = —1) and 
G (0, +oo) for ^ > (s = +1). 

It is worth of noting that the limit (j) — > ±0 (/' — > 1) corresponds to 
the transition to a linear theory: f{R) — > R — 2^^ and R ^ R. This is 
general feature of all nonlinear models f{R)- For example, in the case (13.91 ) 
one obtains f{R) = R{2- exp{A^)) - 2Ad — > R - 2Ad for 0^0. From 
other hand, for particular model (13.91 ). eq. (13.101 ) shows that the point = 
maps into infinity R,R = ±oo. Thus, in this sense, one shall refer to the 
point = as singularity. 

For the model (13.91 ). potential (13.31 ) U{(f)) reads 

^(0) = ^e-^^ (2qs^/\j7\Vse^<f>-s + 2A^) . (3.11) 

It is well known (see e.g. |351ISSII5^) that in order to ensure a stabilization 
and asymptotical freezing of the internal space Mi, the effective potential 
(13.71 ) should have a minimum with respect to both scalar fields ip and (j). Let 
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US remind that the minimum position is chosen with respect to ip ai (f = 0. 
Additionally, the eigenvalues of the mass matrix of the coupled {ip, 0)-field 
system, i.e. the Hessian of the effective potential at the minimum position. 



J := 



(3.12) 



extr 



should be positive definite (this condition ensures the positiveness of the mass 
squared of scalar field excitations). According to the Silvester criterion this 
is equivalent to the condition: 



Ju > , J22 > , det(J) > 



(3.13) 



It is convenient in further consideration to introduce the following nota- 
tions: 

00 := (t)\extr , X := Vse^<^o - s > ^ < 1 . (3.14) 

Then potentials ?7(0), Ueff{f, 4>) and derivatives of the Ueff at an extremum 
(possible minimum) position [ip = 0, 0o) can be rewritten as follows: 

(3.15) 
(3.16) 



t/o = U\e,tr = (1 + sX^) (qs^/\lt]X + A^) 



1 



Uefflextr = "^^l + M^) + h {l + sX^) 



2\c/A 



a + 6 



f^efflextr 



Ri + bUoiX) 



+ {adi + b)h{l + sXy^^ = 



(3.17) 



d^Uefflextr = ch {l ^ sX^'^ - BUq{X) 

= 0, 



2X 



(3.18) 



^(fiifi^ef f\extr 



[a + bf 



Ri + b^UoiX) 



+ {adi + b)'^h{l + sX^) 



(3.19) 
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dl^UeffUr = chadi (1 + sX^f^ , (3.20) 

dl^UeffUr = ch{c-A + 2B){l + sXy^^ 
+ B{A-B)Uo{X) 



{^ + sXy-^'\ (3.21) 

The most natural strategy for extracting detailed information about the lo- 
cation of stability region of parameters in which compactification is possible 
would consist in solving (13.181 ) for X with subsequent back-substitution of the 
found roots into the inequalities (13.131 ) and the equation (13.171 ). To get the 
main features of the model under consideration, it is sufficient to investigate 
two particular nontrivial situations. Both of these cases are easy to handle 
analytically. 

3.3. Case Ag// = 



It can be easily seen from eqs. (13.161 ) and (13.171 ) that condition Ae// = 
Ueff\extr = rcsults in relations 

Ri = 2dih (1 + sX^y^"^ = -^^Uo(X) , di > 2 , (3.22) 

di — 1 



which enables to get from eq. (13.181 ) quadratic equation for X 



{di + 1)^2 + qsdizX - s{di - 1) = , z = 2A^/ ylH (3-23) 
with the following solutions: 



di I I _ / 2 I /I '^i 



= ^^2(dm)p+^V " 

p = ±1. (3.24) 

In the case s = — 1 parameter z should satisfy condition \z\ > zq = 
2 — l/di < 2 and for z = zq two solutions Xp degenerate into one: 
Xp = Xo = -qs^idi-l)/idi + l). 
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Because of conditions h > and e^^° = 1 + sX"^ > 0, the relations 
( 13.221 ) show that parameters Ri and Uo{X) should be non-negative: Ri > 
0, Uq{X) > 0. Obviously, only one of the solutions ( 13.241 ) corresponds to a 
minimum of the effective potential. With respect to this solution let us define 
parameters in the relation ( 13.221 ). Therefore, one musts distinguish now which 
of Xp corresponds to the minimum of U^ff. Let us investigate solutions ( 13.241 ) 
for the purpose of their satisfactions to conditions e^"^" > 0, Uo{X) > and 
Xp > 0. 

The condition e^*^" = 1 + sX^ > 0: 

Simple analysis shows that solutions Xp satisfy this inequality for the 
following combinations of parameters: 



^ > 

(s = +l) 


J9 = ±1 : 2;G(— oo,+oo) 


li<0 
{s = -l) 


j9=+l : 2 G (-2, -zo) U (20, +oo) 
p=-l : z e (-00, -zo) U (zo, 2) 

Z = Zq 



(3.25) 



The condition Uo{X) > 0: 

As appears from eq. ( 13.151 ). this condition takes place if Xp satisfies in- 
equality 2qsXp + z > which leads to the conditions: 



/i > 
(s = +1) 


j9 = +1 : 2G(— oo,+oo) 


/i < 

{s = -l) 


p = +1 : z e {zq, +00) 
p= -1 : z e (zo, 2] 

Z = Zq 



(3.26) 



The condition X„ > 0: 
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This condition is satisfied for the combinations: 



(5 = +1) 


q = +l: < 
q = -1 : < 


z<0: X+>0 X <{) 
z>0: X+>0 X_<0 
\<{): X+<0 X_>0 
z>0: X+<0 X_>0 


li< 
{s = -l) 


g = +1 : < 
q = -l:< 


z<0: X+<0 X_<0 
z>0: X+>0 X_>0 
'z<0: X+>0 X_>0 
z>0: X+<0 X_<0 



(3.27) 



The comparison of (13.251 ). (13.261 ) and (13.271 ) shows that they are simulta- 
neously satisfied only for the following combinations: 



11 > Q 
{s = +l) 


J9 = +1 : g = +1 : 2G(— oo,+oo) 


{s=-l) 


J9=+1 : g = +1 : zG(zo,+oo) 
p=-l : q = +l : z ^ {zq,2) 
q = +1 : z = Zq 



(3.28) 



Additionally, the extremum solutions Xp should correspond to the mini- 
mum of Ueff. The inequalities (13.131 ) are the sufficient and necessary condi- 
tions for that. Let us analyze them in the case of four-dimensional external 
space Dq = 4. Taking into account definitions (13.41 ). (13.81 ). (I3.19I )-( IX^ and 
relations (13.221 ). for Jn, J22 and J21 one gets respectively: 



di + 2 



(3.29) 



J22 — 



6di{2qsXp^z) (l + sX2)" 
{di-l){di + 2){di + 3) 



' rfl+2 



+ 



qs{di + 2) 



(3.30) 



4X3 (di + 3) 



(1 + sX'p) 



+2 



J21 = 



X 



{di-4:){2qsXp + z] 
(di- l)(di + 2) 

(l + sX^)-^ 



+2 



di{di + 3) 



(3.3i; 
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It is supposed in these equations that each of Xp can define zero minimum 
of Ueff. In what follows, one shall check this assumption for every Xp with 
corresponding combinations of signs of the parameters s and q in accordance 
with the table (ESHD- 

According to the Silvester criterion ( 13.131 ). Jn should be positive. Thus 
eqs. (13.221 ) and (13.291 ) result in the following conclusions: the potential Uq 
should be positive Uq > 0, the internal space should have positive curvature 
Ri > (hence, di > 1) and its stabilization (with zero minimum Ag// = 0) 
takes place only in the present of form-field {h > 0). Transition from the 
non-negativity condition Uq > to the positivity one Uq > corresponds to 
the only substitution in (13.261 ) (zq, 2] — > (^o, 2) for the case s = —1, p = — 1. 
Exactly this interval {zo,2) appears in concluding table (13.281 ). Therefore, 
Jii is positive for all Xp from the table (13.281 ). 

Concerning expressions J22 and det( J) = Jn J22 — graphical plotting 
(see Fig J3^ and Fig J3.5l ) demonstrates that they are negative for s = +1, p = 
+1, q = +1 and s = —1, p = —1, q = +1 but positive in the case s = 
— 1, p = +1, q = +1. For this latter combination z G {zq^oo). The case 
s = — 1, q = +1 and z = zq should be investigated separately. Here, Xp = 
Xq = {d\ — 1)/ {d\ + 1) and for J 22 and J21 one obtains: 

f 12(d,-l)-(d, + 2)^ \ 

V(rfi + i)2(rfi + 2)(di + 3);' ^ ^ ^ 

+2 



J2x= - (1 + 

^ / 2v/2(d,-4)(rf,-i; 

V(di + l)2(di + 2)v'^ii(c/i + 3) ' 

It can be easily seen from eqs. (13.321 ) and (13.331 ) that J22 > for c/i 7^ 4 and 
J22 = J21 = for o?i = 4. Additionally, det( J) > for dx ^ 4. 

Thus, it can finally concluded that zero minimum of the effective potential 
Ueff takes place either for s = —1, q = +1, z G (zq, 00), Vdi > 1 (position of 
this minimum is defined by solution (13.241 ) with p = +1) or for s = — 1, q = 
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+1, z = zq,\/ di 7^ 4. Concerning the signs of parameters, one obtains that 
/i < and Ad > 0. 




Figure 3.4: Typical form of J22/ ^/\JA (^Q- (I3.30p ) for parameters s = +1, q = +1, p = +1 
(left) and s = —1, q = +1, p = —1 (right). For both of these combinations of the parameters, 
J22 < 0. 



J22/V17r s = -l p= + l q= + l Det(J)/|^i| s._i p. + i q=+i 




zo(di) Zo(di) 



Figure 3.5: Typical form of J22/ ^AiA (left) and det{J)/\fi\ (right) (see eqs. fl3:29|) - fl3T[]) ) 
for parameters s = —1, q = +1, p = +1. For this combination of the parameters, J22 > 
and det J > 0. 



3.4. Decoupling of excitations: di = Dq 

It can be easily seen from eq. (13.81 ) that in the case di = Dq parameter 
c = that leads to condition d'^^Uef f\extr = (see eq. (13.201 ). Thus, the 
Hessian (13.121 ) is diagonalized. It means that the excitations of the fields ip 
and (j) near the extremum position are decoupled from each other 3. 



'''In the vicinity of a minimum of the effective potential, squared masses of these excitations are = Jn 
and ml = J22. 
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Dropping the h term in eq. (13.181 ) (because of c = 0) and taking into 
account eq. (13.151 ). one obtains quadratic equation for X 



[D + 2)X2 + DqszX - (L» - 2)s = , 



(3.34) 



which for Dq = di exactly coincides with eq. (12.431 ). Thus, in spite of the 
fact that the condition Ag// = is not apphed directly, one obtains in the 
case di = Dq precisely the same solutions (13.241 ). However, parameters Uq 
and h satisfy now relations different from (13.221 ). For example, for the most 
physically interesting case Dq = di = 4, eqs. (13.161 ) and (13.171 ) result in the 
following relations: 



Ri 



\u,{x) + h 



h 



Nonzero components of the Hessian read 



Jn = ^[9h 



J22 



Uq{X)] , 

4(z + 2gsXp) ^ 3gs(l + sX2)2/3' 



(3.35) 



(3.36) 



14X3 



21(l + sX|)4/3 

Now, let us to look for a positive minimum of the effective potential. It 
means that Ag// > 0, Jvi > and J22 > 0. From the positivity of Jn and 
Ag// follows respectivelyl^: 



Jn > : 16/i > Ri> IWq{X)/^ > 8A, 



eff 



(3.37) 



and 



Aeff > : h> Ri/16 > Uq{X)/9 > h/3>0 . (3.38) 

These inequalities show that for the considered model positive minimum of 
the effective potential is possible only in the case of positive curvature of the 
internal space Ri > and in the presence of the form field {h > 0). 

To realize which combination of parameters s^p and q ensures the mini- 
mum of the effective potential, one should perform analysis as in the previous 



di — 4) relations p.34p . Jn in (|3.36l) and inequalities 
coincide with the analogous expressions in paper |58j with quadratic nonlinear model. This is 



^It is interesting to note that (in the case Do 

(ESZD, (E 



not surprising because they do not depend on the form of nonlinearity f{R) (and, consequently, on the form 
of U[4>)). However, the expressions for J22 are different because here the exact form of J7(0) is used. 
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case with Ag// = 0. However, there is no need to perform such analysis here 
because solutions of eq. (13.341) coincides with (13.241) and all conditions for 
Xp and Uq are the same as in the previous section. Thus, one obtains con- 
cluding table of the form (13.281 ). Additionally, it can be easily seen that 
expressions of J22 in (13.361) and (13.301 ) exactly coincide with each other if one 
put = 4 in the latter equatioijj. Hence, Fig J3.4l and Fig J3.5l can be used 
(for the lines with di = 4) for analyzing the sign of J22- With the help of 
these pictures as well as keeping in the mind that J22{d\ = 4, z = zq) = 0, 
follows that the only combination which ensures the positive minimum of 
f/g// is: s = — 1 , j3 = +1 , q = +1 and z ^ {zq = 1.936, +00). It is clear that 
potential Uq{X) in eqs. (I3.35I) -( IX^ is defined by solution of eq. (13.341 ) (i.e. 
(13.241) for di = 4) with this combination of the parameters. Because s = — 1 
and 2 > , the parameters /i and Kj:, should have the following signs: /i < 
and > 0. 

Additionally, it is easy to verify that second solution of eq. (13.241) X 
(with p = —1 and s = — 1,^ = +!) does not correspond to the maximum 
of the effective potential Ueff. Indeed, in this case d(j)Ueff{X ) = but 
d^U,ff{X_) ^ 0. 

Fig. 13.61 demonstrates the typical profile = of the effective potential 
Ueff{<-P: (j)) in the case of positive minimum of Ueff considered in the present 
section. This picture is in good concordance with the table (13.281) . According 
to this table, positive extrema of t/g// are possible only for the branch q = +1 
of the solution (13.101) (solid lines in Fig. 13.61 ). It is obvious that for < 2; < 2 
one can have 3 such extrema: one for positive /i>0— >s = +l and two for 
negative fi < ^ s = —1. These investigations show that in the left half 
plane (i.e. for u < 0) the right extremum [p = +1 in eq. (13.241) ) is the 
local minimuno and the left maximum {p = —1) is not the extremum of 
Ueff because here d^Ueff 7^ 0. Analogously, maximum in the right half plane 



^It follows from the fact that in (I3.30p Dq — 4. Although in this equation the relation (I3.22p between h 
and Uq is used, it enters here in the combination which is proportional to c. Thus, this combination does 
not contribute if additionally c?i — 4. 

^°It can be easily seen for the branch q = +1, s = —1 that {/{(p — >■ — oo) — ?> +oo for z > 2 and 
U{(f> —7- — oo) — > — oo for z < 2. Thus, for z > 2 this minimum becomes global one. 
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/i > (which corresponds to {p = +l)-solution ( 13.241 )) is not the extremum of 
Ueff. For completeness of picture, hnes corresponding to the branch q = —1 
(dashed hnes in Fig. 13.61 ) are also included. The minimum of the right dashed 
line (for Xp with s = +l, j9=— 1, q = —1) does not describe the extremum 
of Ueff because again d^pUeff ^ 0. 



Ueff 
1.^ 




Figure 3.6: Profile = of the effective potential Ueff{^, 0) for parameters z = 1.99 (A^ = 
1.99/2 , \fi\ = 1), Dq = di = 4 and h = Uq/6. The rest of parameters can be found from 
relations (I3.35p . Solid and dashed lines describe branches q = +1 and q = —1 respectively. 



For the solid line, there is only one local minimum of Ueff which is defined by solution Xp 
with the following parameters: s = —1, p = +1, q = +1. Left {s = —1, p = —1, q = +1) 
and right {s = +1, p = +1, q = +1) maxima are not extrema of Ueff because in these 
points d^Ueff = but d^Ueff 7^ 0. Analogously, for right dashed line the solution Xp (with 
s = +l, p = —1, q = —1) does not correspond to the extremum of Ueff- 



3.5. Cosmic acceleration and domain walls 

Let us consider again the model with di = Dq = 4 in order to define 
stages of the accelerating expansion of our Universe. It was proven that for 
certain conditions (see (I3.35I )-( IX^ ) the effective potential Ueff has local (for 
zo < 2; < 2) or global (for z > 2) positive minimum. The position of this 
minimum is ((/9 = 0, = (1/^) ln(l + sX^)) where s = —1 , p = +1 , g = +1 
and z G {zq^+oo). Obviously, positive minimum of the effective potential 
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plays the role of the positive cosmological constant. Therefore, the Universe 
undergoes the accelerating expansion in this position. Thus, one can "kill two 
birds with one stone" : to achieve the stable compactification of the internal 
space and to get the accelerating expansion of our external space. 

Let us associate this acceleration with the late-time accelerating expansion 
of our Universe. As it follows from eqs. (13.351 ) and ( 13.381 ). positive minimum 
takes place if the parameters are positive and the same order of magnitude: 
Ae// ~ i?i ~ U{X) ^ h > 0. On the other hand, in KK models the size of 
extra dimensions at present time should be 6(o)i ^ lO^^^cm ~ ITeV^^. In 
this case Ri > b^^^^ - lO^^cm-^ Thus, for the TeV scale of 6(o)i ~ ITeV 
one gets that Ag// ~ i?i ~ U{X) ~ /i ~ ITeV^. Moreover, in the case of 
natural condition A^) ~ \/JiA it follows that the masses of excitations ~ 
rricj) ~ ITeV. The above estimates clearly demonstrate the typical problem of 
the stable compactification in multidimensional cosmological models because 
for the effective cosmological constant one obtains a value which is in many 
orders of magnitude greater than observable at the present time dark energy 
~ 10~^^cm~^. The necessary small value of the effective cosmological constant 
can be achieved only if the parameters Ri , U {X) , h are extremely fine tuned 
with each other to provide the observed small value from equation Ae//(X) = 
Uo{X)/3 — h. Two possibilities to avoid this problem can be seen. Firstly, 
the inclusion of different form- fields/ fluxes may result in a big number of 
minima (landscape) [HllISS] with sufficient large probability to find oneself 
in a dark energy minimum. Secondly, the restriction Ri ~ 6^^-^ ~ lO^^cm^^ 
can be avoided if the internal space is Ricci-flat: Ri = 0. For example, the 
internal factor-space Mi can be an orbifold with branes in fixed points (see 
corresponding discussion in [86]). 

The WMAP three year data as well as CMB data are consistent with wide 
range of possible inflationary models (see e.g. |[I3])- Therefore, it is of interest 
to get the stage of early inflation in the given model. It is well known that 
it is rather difficult to construct inflationary models from multidimensional 
cosmological models and string theories. The main reason of it consists in the 
form of the effective potential which is a combination of exponential functions 
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(see e.g. eq. (13.71 )). Usually, degrees of these exponents are too large to result 
in sufficiently small slow-roll parameters (see e.g. |I3Z])- Nevertheless, there 
is a possibility that in the vicinity of maximum or saddle points the effective 
potential is flat enough to produce the topological inflation |H7H93] . Let us 
investigate this possibility for the given model. 

As stated above, the value cp = corresponds to the internal space value 
at the present time. Following this statement, the minimum of the effective 
potential at this value of (p is found. Obviously, the effective potential can 
also have extrema at ip 0. Let us investigate this possibility for the model 
with di = Dq = 4, i.e. for c = 0. In this case, the extremum condition of the 
effective potential reads 

^.^e//|,„,^„ = -^i?i(a + 6)e(«+^)^o + 6[/oe^^« 

+ (a(ii + 6)/ie("^i+^)^" = (3.39) 



and 



dU 

= 0; =^ 



= . (3.40) 



Here, <pq and 0o define the extremum position. It clearly follows from eq. 
(13.401 ) that 00 is defined by equation dU/dcj) = which does not depend 
on ip. Therefore, extrema of the effective potential may take place only for 
00 which correspond to the solutions Xp of eq. (13.341 ) and different possible 
extrema should lie on the sections Xp = const. So, let us take X+ (with 
s = — 1 and g = +1) which defines the minimum of Ueff in the previous 
section. Hence, Uq in eq. (13.391 ) is the same as for eq. (13.351 ). 

Let us define now (/9o from eq. (13.391 ). With the help of inequahties (13.371 ) 



and (13.381 ) one can write h = hUq where n G (1/9,1/3). Taking also into 
account relations (13.351 ). eq. (13.391 ) can be written as 



where the definition is done y = exp (a(/?o) = exp ((^o/VS) and put di = Dq = 
4. Because y = 1 \s the solution of eq. (13.411 ). the remaining three solutions 



50 



satisfy the following cubic equation: 



y^ + y^ + y 



1 _^ 

3n 



(3.42) 



It can be easily verified that the only real solution of this equation is 



2v + - 

V 



(3.43) 



where 



21/3 



n 



V = 



(9n2 + 7n3 + 3V9n4 + 14^6 + 9^6) 
1 1' 
9' 3 



1/3 



(3.44) 



Thus (po{yo) and 4)o{X+) define new extremum of U^ff. To clarify the type of 
this extremum one should check signs of the second derivatives of the effective 
potential in this point. First of all let us remind that in the case c = mixed 
second derivative disappears. Concerning second derivative with respect to 
0, one obtains 



J22 — 



= e 



because in previous section was obtained d'^U /dcj)'' 
five with respect to (p reads 



> (3.45) 
> 0. Second deriva- 



Jn = 



eff 



ylUo 



-6 Q + n) + J 



+ Unyi] = 2ylUo 



{3n+l)yo - - 



(3.46) 



where eq. (13.411 ) is taken into account. Simple analysis shows that (3n + 
l)yQ < 4/3 for n G (1/9, 1/3). Keeping in mind that Uq > one obtains Jn < 
0. Therefore, the extremum is the saddle surfac^. Figure [231 demonstrates 
contour plot of the effective potential in the vicinity of the local minimum 
and the saddle point. 



Similar analysis performed for the branch with s = +1, q = —1 (right dashed line in the Fig. I3.6P shows 
the existence of the global negative minimum with f =/= along the section X_ = const. 
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Figure 3.7: Contour plot of the effective potential Ueff{^, 0) for parameters z = 1.99 {Ad = 
1.99/2 , = 1), Dq = di = 4 and h = Uq/G. The rest of parameters follows from relations 
(13.351) . The branch corresponding tos = — l,g = -|-lis chosen. This plot clearly shows the 
minimum and the saddle points of the effective potential. 

Therefore, very interesting possibility become possible for the production 
of an inflating domain wall in the vicinity of the saddle point. The mechanism 
for the production of the domain walls is the following [HZIIHB]- If the scalar 
field (p is randomly distributed, some part of the Universe will roll down to 
(p = 0, while in others parts it will run away to infinity. Between any two such 
regions there will appear domain walls. In Ref. [i89], it was shown for the case 
of a double- well potential Vdw{^) = (-^/4)(</?^ — ^^)^ that a domain wall will 
undergo inflation if the distance ^ between the minimum and the maximum 
of Vdw exceeds a critical value ^cr = O.SSMpi — > Ko^cr = 1-65. In this case 
it means that the distance \(fQ\ between the local minimum and the saddle 
point should be greater than ^cr- |</^o| ^ 1-65. Unfortunately, for this model 
(po{n) < 1.65 if n G (1/9, 1.3). For example, in the most interesting case n — > 
1/3 (where Ae// — > (see eq. (13.351 ))) one obtains \ipo\ — > 1.055 which is less 
than ^cr- Moreover, the domain wall is not thick enough in comparison with 
the Hubble radius. The ration of the characteristic thickness of the wall to the 
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horizon scale is given by Vy^H « |?7e///35^^?7e//|^^ (j>{x+) ~^ 0.454 for n — > 1/3 
which is less than the critical value 0.48 for a double- well potential. Thus, here 
there is no a sufficiently large (for inflation) quasi-homogeneous region of the 
energy density. And the given potential is too steep. Obviously, the slaw roll 
parameter e « (1/2) (dipUeff/Ueff)^^^^ is equal to zero in the saddle point. 
However, another slow roll parameter t] « \d'^^Ueff/Ueff\^ ^ —1.617 for 
n — > 1/3. Therefore, the domain walls do not inflate in contrast to the case 
in Ref. PDHSH. 

In Fig. 13.81 the comparison is presented between the potential (solid line) 
and a double-well potential (dashed hue) in the case n = {1 — 0.001)/3. 
It is shown that the given potential is flatter than a double-well potential 
around the saddle point. However, calculations show that it is not enough 
for inflation. 



Ueff (V,(Po) 




Figure 3.8: Comparison of the potential Ueff{f, 0o) with a double well potential for param- 
eters z = 1.99 {Ad = 1.99/2 , = 1), Dq = di = A and n = {1 - 0.001)/3. 

Summary 

It was shown that positive minimum of the effective potential plays the 
double role in this model. Firstly, it provides the freezing stabilization of 
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the internal spaces which enables to avoid the problem of the fundamen- 
tal constant variation in multidimensional models [ESIISS]- Secondly, it en- 
sures the stage of the cosmic acceleration. However, to get the present-day 
accelerating expansion, the parameters of the model should be fine tuned. 
Maybe, this problem can be resolved with the help of the idea of landscape 
of vacua [B11IH5]. 

It was additionally found that given effective potential has the saddle 
point. It results in domain walls which separates regions with different vacua 
in the Universe. These domain walls do not undergo inflation because the 
effective potential is not flat enough around the saddle point. 

It is worth of noting that minimum in Fig. 13.61 (left solid line) is metastable. 
In other words, classically it is stable but there is a possibility for quantum 
tunnelling both in and in if directions (see Fig. 13.71 ). One can avoid this 
problem in direction in the case of parameters z > 2 (see footnote [TUj ). 
However, tunnelling in (p direction (through the saddle) is still valid because 
Ueff{(f,(l)o) ~ e^'^U{(J)o) — > for — > — oo which is less than any positive 
Ae//. It may result in the materialization of bubbles of the new phase in the 
metastable one (see e.g. [ST]). Thus, late-time acceleration is possible only 
if characteristic lifetime of the metastable stage is greater than the age of 
the Universe. Careful investigation of this problem (including gravitational 
effects) is rather laborious task which needs a separate consideration. As it 
was mentioned in footnote [lU there is also the global negative minimum for 
right dashed line in Fig. 13.61 (it corresponds to the point {i^ = 0.67, = 
1.66) for parameters taken in Fig. 13.71 ). This minimum is stable both in 
classical and quantum limits. However, the acceleration is absent because of 
its negativity. 

Another very interesting feature of the model under consideration consists 
in multi-valued form of the effective potential. As it can be easily seen from 
eqs. (13.71 ) and (13.111 ). for each choice of parameter /i potential U{(J)) (and 
consequently Ue/f) has two branches {q = ±1) which joint smoothly with 
each other at = (see Fig. 13.61 ). It gives very interesting possibility 
to investigate transitions from one branch to another one by analogy with 
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catastrophe theory or similar to the phase transitions in statistical theory. 
However, as it was mentioned above, in this particular model the point = 
corresponds to the singularity R ^ ±00. Thus, the analog of the second 
order smooth phase transition through the point = is impossible in this 
model. Nevertheless, there is still a possibility for the analog of the first order 
transition via quantum jumps from one branch to another one. 

To complete, let us investigate some limiting cases. Firstly, the limit /i — > 
(for arbitrary Dq and di) where the form-fields are absent. From eqs. (13.151 ) 
- (13.211 ) one obtains the following system of equations: 

fi, = JLc7„(X). U.„U^^^U„iX) (3.47) 

and 



J,, = -abUoiX), J2i = 0, 



J22 — 



B{A-B)-'^{l + eX'- 



U,>(X) . (3.48) 



Since for minimum should hold true the condition Jn > 0, one arrives at 
the conclusion: Uq^ Uejf < 0. Consequently, the minimum of the effec- 
tive potential as well as the effective cosmological constant is negative and 
accelerating expansion is absent in this limit. Therefore, the presence of the 
form-fields is the necessary condition for the acceleration of the Universe in 
the position of the freezing stabilization of the internal spaces. Additionally, 
it can be easily seen that the extremum position equation takes the same form 
as (13.341 ). Simple analysis show that minimum takes place for the branch: 



s = +1 (i.e. fi > 0)^ p = —1, q = —1 and z G (—00, +00). If additionally is 
imposed that 2; — > (i.e. Ajj ^ and /i is fixed) then the results of Ref. jS^ 
are reproduced. 



55 



CHAPTER 4 

INFLATION IN MULTIDIMENSIONAL 
COSMOLOGICAL MODELS 

In this chapter, a multidimensional cosmological models are considered with 
linear, nonlinear quadratic and quartic actions with a monopole form 
field as a matter source; and pure gravitational model with B? + nonlin- 
earities. The inflation is investigated in these models. It is shown that 
and R^ models can have up to 10 and 22 e-foldings, respectively. These val- 
ues are not sufficient to solve the homogeneity and isotropy problem but big 
enough to explain the recent CMB data. Additionally, model can provide 
conditions for eternal topological inflation. However, the main drawback of 
the given inflationary models consists in a value of spectral index which is 
less than observable now « 1. For example, in the case of R^ model one 
finds Us ~ 0.61. 

For the model R^ + . the effective scalar degree of freedom (scalaron) 
has a multi- valued potential U (0) consisting of a number of branches. These 
branches are fitted with each other in the branching and monotonic points. 
In the case of four-dimensional space-time, it is shown that the monotonic 
points are penetrable for scalaron while in the vicinity of the branching points 
scalaron has the bouncing behavior and cannot cross these points. Moreover, 
there are branching points where scalaron bounces an infinite number of times 
with decreasing amplitude and the Universe asymptotically approaches the 
de Sitter stage. Such accelerating behavior is called bouncing inflation. For 
this accelerating expansion there is no need for original potential U (0) to 
have a minimum or to check the slow-roll conditions. A necessary condition 
for such inflation is the existence of the branching points. This is a new 
type of inflation. It is shown that bouncing inflation takes place both in the 
Einstein and Brans-Dicke frames. 
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4.1. Linear model 



To start with, let us define the topology of the given models and consider 
a factorizable D-dimensional metric 

9<°)(,-.) + L|,e^^'(V), (4.1) 

which is defined on a warped product manifold M = Mq x Mi. Mq describes 
external Do-dimensional space-time (usually Dq = 4) and Mi corresponds to 
di-dimensional internal space which is a flat orbifolco with branes in fixed 
points. Scale factor of the internal space depends on coordinates x of the 
external space-time: ai{x) = Lpie^^^^\ where Lpi is the Planck length. 

First, let us consider the linear model f{R) = R with D-dimensional action 
of the form 



S = 



^ / d^'x^l^l -2An] + + Sb . (4.2; 



A.£) is a bare cosmological constanlo. In the spirit of Universal Extra Di- 
mension models P51 - I102j . the Standard Model fields are not localized on the 
branes but can move in the bulk. The compactification of the extra dimen- 
sions on orbifolds has a number of very interesting and useful properties, e.g. 
breaking (super) symmetry and obtaining chiral fermions in four dimensions 
(see e.g. paper by H.-C. Cheng at al in pg] - ll02 j). The latter property gives a 
possibility to avoid famous no-go theorem of KK models (see e.g. fl03lll04] ). 
Additional arguments in favor of UED models are listed in jl05l . 



Following a generalized Freund- Rubin ansatz jBD] to achieve a spontaneous 
compactification M — )► M = Mq x Mi, one endows the extra dimensions with 
real-valued solitonic form field F^^^ with an action: 

'- = -\L/^^\hM' (4.3) 



12 



For example, 5*^/^2 and T"^ /Z2 which represent circle and square folded onto themselves due to Z2 
symmetry. 

^■^Such cosmological constant can originate from D-dimensional form field which is proportional to the 



Z3-dimensional world-volume: p^^ - Q = {C / \/\cp^\)e'^'^'^ ■■■'-^ . In this case the equations of motion gives 
C = const and term in action is reduced to {1/D\)Fmn...qF^^^---'^ = -C^. 
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This form field is nested in (ii -dimensional factor space Ml, i.e. IS propor- 
tional to the world- volume of the internal space. In this case (l/di!) (i^^^^) = 
fi/a^i^ , where /i is a constant of integration 

Branes in fixed points contribute in action functional ( 14.21 ) in the form 



JMn 



fixed 
points 



fixed 
point 



(4.4) 



where g^^\x) is induced metric (which for the given geometry (14.11) coincides 
with the metric of the external space-time in the Brans-Dicke frame) and 
is the matter Lagrangian on the brane. In what follows, that the case where 
branes are only characterized by their tensions L;,(/j) = ~^k) , k = 1,2,. ..,m 
is considered, where m is the number of branes. 

Let /3q be the internal space scale factor at the present time and = (3^ — 
13^ describes fiuctuations around this value. Then after dimensional reduction 
of the action (14.11) and conformal transformation to the Einstein frame gj^J = 

' giiJ , one arrives at effective Do-dimensional action of the form 




- ~g^'^^''d,ipd,ip-2U,ff{ip)] , (4.5) 

where scalar field ip is defined by the fiuctuations of the internal space scale 
factor: 

and G := kI/Stt := Kjj/{87TVd^) {Vd^ is the internal space volume at the 
present time) denotes the Do-dimensional gravitational constant. The effec- 
tive potential Ueff{(p) reads (hereafter Dq = 4): 



2di 
T2 ^ 



Ae V-ii+s 



2(ii 

— V 



(4.7) 



where /f = i4,fi/alo)i and A = -kI Y,'k=i nk)- 



58 



Now, this potential should be investigated from the point of the external 
space inflation and the internal space stabilization. First, it is clear that 
internal space is stabilized if Ueff{(f) has a minimum with respect to (f. 
The position of minimum should correspond to the present day value (p = 
0. Additionally, it needs to be demanded that the value of the effective 
potential in the minimum position is equal to the present day dark energy 
value Ueff{(f = 0) ^ Ade ^ lO^^^cm^^. However, it results in very flat 
minimum of the effective potential which in fact destabilizes the internal 
space [HS]- To avoid this problem, the case of zero minimum Ueff{ip = 0) = 



shall be considered. 

The extremum condition dUefj/difl^^Q = and zero minimum condition 
Ueff{<-p = 0) = result in a system of equations for parameters A^:;, and A 
which has the following solution: 

A^ = = A/2 . (4.8) 

For the mass of scalar field excitations (gravexcitons/radions) one obtains: 
= (fUeff/d(p^\^^^ = (4di/(di + 2))Az5. In Fig. the effective po- 
tential (14.71) is presented in the case = 3 and Kb = 10. It is worth of 
noting that usually scalar fields in the present paper are dimensionlesqlj and 
?7e//, Ao, fi, X are measured in Mpi units. 

Let us turn now to the problem of the external space inflation. As far as 
the external space corresponds to our Universe, the metric ^(0) is taken in 
the spatially flat Friedmann-Robertson- Walker form with scale factor a{t). 
Scalar field depends also only on the synchronous/cosmic time t (in the 
Einstein frame). 

It can be easily seen that for cp >> (more precisely, for (p > (pmax = 
^y{dl+~2^)J2dl In 3) the potential (14.71) behaves as 

Ueffip) ^ Adc-^^"^ ^ (4.9) 

with 

2di 



*To restore dimension of scalar fields one should multiply their dimensionless values by Mpi/\/8tt. 
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Figure 4.9: The form of the effective potential ( 14. 7p in the case di = 3 and Ad = /f = A/2 
10. 



It is well known (see e.g. fl06H109j ) that for such exponential potential scale 
factor has the following asymptotic form: 



a{t) - t^/g 



(4.11) 



Thus, the Universe undergoes the power-law inflation if g < 2. Precisely this 
condition holds for eq. (14.101 ) if di > 1. 

It can be easily verified that (p > (pmax is the only region of the effec- 
tive potential where inflation takes place. Indeed, in the region (p < the 
leading exponents are too large, i.e. the potential is too steep. The lo- 
cal maximum of the effective potential Ueff\^^^ = (4/27) A^) at (pmax = 
a/ {di + 2)/2di In 3 is also too steep for inflation because the slow-roll pa- 



rameter r]r^ 



_ 1 d^Ueff 
Ueff dip^ 



_ Ml 



di+2 



1 < \ilmax\ < 3 and does 
Topological inflation is also 



not satisfy the inflation condition \r]\ < 1 
absent here because the distance between global minimum and local max- 
imum (pmax = \/ {di + 2)/2di In 3 < 1.35 is less than critical value (pcr > 1.65 
(see [SlllSSllHS])- It is worth of noting that ijmax and cpmax depend only on 
the number of dimensions di of the internal space and do not depend on the 
hight of the local maximum (which is proportional to Ajj). 

Therefore, there are two distinctive regions in this model. In the first 
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region, at the left of the maximum in the vicinity of the minimum, scalar 
field undergoes the damped oscillations. These oscillations have the form of 
massive scalar fields in our Universe (in jllOj these excitations were called 
gravitational excitons and later (see e.g. jlllj ) these geometrical moduli os- 
cillations were also named radions). Their life-time with respect to the decay 
— > 27 into radiation is p5llll2l[TT5] r ~ {Mpi/m)^Tpi. For example, in 
given case one obtains r ~ 10 s, 10^^ s for m ~ 10 TeV, 10^ TeV correspond- 
ingly, where = (4o?i/ {di+2))Ajj. Therefore, this is the graceful exit region. 
Here, the internal space scale factor, after decay its oscillations into radia- 
tion, is stabilized at the present day value and the effective potential vanishes 
due to zero minimum. In second region, at the right of the maximum of the 
potential, our Universe undergoes the power-low inflation. However, it is 
impossible to transit from the region of inflation to the graceful exit region 
because given inflationary solution satisfies the following condition (p > 0. 
There is also serious additional problem connected with obtained inflation- 
ary solution. The point is that for the exponential potential of the form (14.91 ). 
the spectral index reads jl06L lT08] 



15 



n. = |-^. (4.12) 

z q 

In the case (14.101 ). it results in Us = I — di. Obviously, for di > 1 this value is 
very far from observable data ~ 1. Therefore, it is necessary to generalize 
given linear model. 



4.2. Quadratic model 



One of possible ways for generalizing of the effective potential making it 
more complicated and having more reach structure is an introduction of an 
additional minimal scalar field 0. It is possible to do "by hand" , inserting 

^^With respect to conformal time, solution (j4.11l) reads a(rj) ^ ry^+^ where /3 — —(4 — q)/{2 — q). It was 
shown in [114] that for such inflationary solution (with q < 2) the spectral index of density perturbation is 
given hy Us = 2f3 + 5 resulting again in (j4.12p . 
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minimal scalar field (f) with a potential U (0) in linear action ( I4.2l r1. Then, 
effective potential takes the form 



2di 

Ueff{ip,(l)) = e"V^^ 



Ae 



2d^ 



(4.13) 



where = in ( lO ). 

However, it is well known that scalar field can naturally originate from 
the nonlinearity of higher-dimensional models where the Hilbert-Einstein lin- 
ear lagrangian R is replaced by nonlinear one f{R)- These nonlinear theories 
are equivalent to the linear ones with a minimal scalar field (which repre- 
sents additional degree of freedom of the original nonlinear theory). It is not 
difficult to verify (see e.g. [SIIEB]) that nonlinear model 



- Y, [ d'x^lffC'WI nk) (4.14) 



is equivalent to a linear one with conformally related metric 

gio) = e-^*/<--)5ir (4.15) 
plus minimal scalar field (f) = ln[df /dR]/A with a potential 



/7(0) = ie-^^ [i?(0)e^^ - /(i?(0))] , (4.16) 



where 



A = 




{di + 2) 



1) V ^^1 + 3 



5 = , = . (4.17) 

V{D-2){D-1) (di + 2) 



"'^^If such scalar field is the only matter field in these models, it is known (see e.g. [571lll5j ) that the effective 
potential can has only negative minimum, i.e. the models are asymptotical AdS. To uplift this minimum to 
nonnegative values, it is necessary to add form- fields [55] . 
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After dimensional reduction of this linear model, one obtains an effective 
Do-dimensional action of the form 

1 



2^0 J Mo 



g 



(4.18) 



with effective potential exactly of the form ( 14.131 ). It is worth to note that 
it is supposed that matter fields are coupled to the metric g^^^ of the linear 
theory (see also analogous approach in jll6j ). Because in all considered below 
models both fields (p and (J) are stabilized in the minimum of the effective 
potential, such convention results in a simple redefinition/rescaling of the 
matter fields and effective four-dimensional fundamental constants. After 
such stabilization, the Einstein and Brans-Dicke frames are equivalent each 
other (metrics g''^^ and g^^^ coincide with each other), and linear (/*^^) and 



nonlinear g^^^ metrics in (I4.15f 



17 



are related via constant prefactor (models 



became asymptotically linear) 

Let us consider first the quadratic theory 

f{R) = R + ^R^ -2Ad. 
For this model the scalar field potential (14.161 ) reads: 



40 



1) +2Ad 



(4.19) 



(4.20) 



It was proven [1 115] that the internal space is stabilized if the effective 
potential (14.131 ) has a minimum with respect to both fields (p and (J). It can be 
easily seen from the form of Ueff{(p, 4>) that minimum 0o of the potential U (0) 
coincides with the minimum of Ueff{<-p, (j)) : dU/d(j)\^^ = — > d(j)Ueff\^^ = 0. 
For minimum ?7(0o) one obtains |57j : 



-D 



U{^o) = [{xo - 1)2 + 8^A^] , 



{4.2i; 



where the denotation xq := exp{A(f)Q) is applied, taking into account that 
exp(^0o) = - 5 + v^^2 ^ (^2A - B)8sB ^Ad^ /{2A-B). It is the 

"'^^However, small quantum fluctuations around the minimum of the effective potential distinguish these 
metrics. 
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global minimum and the only extremum of U{(j)). Nonnegative minimum 
of the effective potential U^jf takes place for positive ^, Kb > 0. If ^, > 0, 
the potential U{(J)) has asymptotic behavior U{(j)) — > +oo for — > ±oo. 

The relations (14.81 ). where one should make the substitution Ad — > U{(f)o), 
are the necessary and sufficient conditions of the zero minimum of the effective 
potential Ueff{<-p, (j)) at the point {cp = 0,(f) = 0o)- Thus, if parameters of the 
quadratic models satisfy the conditions U{(j)o) = /f = A/2, one arrives at 
zero global minimum: t/e//(O,0o) = 0. 

It is clear that profile = 0o of the effective potential t/g// has a local 
maximum in the region cp > because Ueff{(p^ (J) = 0o) — > if — )• +oo. Such 
profile has the form shown in Fig. 14.91 Thus, the effective potential t/g// has 
a saddle point {(p = (pmax^ = 0o) where (pmax = \/ {di + 2)/2(ii In 3. At this 
point ^e//|„,aa; ~ (4/27)t/(0o)- The Figure STIU] demonstrates the typical 
contour plot of the effective potential (14.131 ) with the potential U (0) of the 
form (14.201 ) in the vicinity of the global minimum and the saddle point. 




Figure 4.10: Contour plot of the effective potential Ueff{^,(f)) fl4.13p with potential f/(0) of 
the form fl4.20p for parameters di = 1 ,^A£) = 1 and relations t/(0o) = /i = '^/2 • This plot 
clearly shows the global minimum and the saddle. The colored lines describe trajectories for 
scalar fields starting at different initial conditions. 
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Let us discuss now a possibility of the external space inflation in this 
model. It can be easily realized that for all models of the form (14.131) in the 
case of local zero minimum at {(f = 0,00), the effective potential will also 
have a saddle point at {ip = i^max-, 0o) with (pmax = \/ {di + 2)/2di In 3 < 
(fcr = 1-65 and the slow-roll parameter |?7(^| in this point cannot be less than 
1: |?7<^| = 3(ii/((ii + 2) > 1. Therefore, such saddles are too steep (in the 
section = 0o) for the slow-roll and topological inflations. However, as one 
shall see below, a short period of de-Sitter-like inflation is possible if one start 
not precisely at the saddle point but flrst move in the vicinity of the saddle 
along the line <^ ~ '^max with subsequent turn into zero minimum along the 
line ~ 00- Similar situation happens for trajectories from different regions 
of the effective potential which can reach this saddle and spend here a some 
time (moving along the line (f ~ ^max)- 

Let us consider now regions where the following conditions take place: 

r, / 2di I 2di 

/7(0)>/iW^'", Ae^V^^. (4.22) 

For the potential (14.201 ) these regions exist both for negative and positive 0. 
In the case of positive with exp(740) ;> max {l , (8^Ai:»)^/^} one obtains 

/7e//«^e-^^^e^^ (4.23) 

where q is deflned by Eq. mUI ). qi := {2A - Bf = dl/Mi + 2){di + 3)] and 
q > qi- For potential (14.231) the slow-roll parameters arellj: 

e « r/i « r/2 « I + |- (4.24) 

and satisfy the slow-roll conditions 6,771,772 < 1. As far as it is known, there 
are no analytic solutions for such two-scalar-fleld potential. Anyway, from 
the form of the potential (14.231) and condition q > qi one can get an estima- 
tion a « with s > 2/q (e.g. 2/q = 3,2,5/3 for di = 1,2,3, respectively). 

-'^^In the case of n scalar fields ipi {i — 1, . . . ,n) with a flat (cr— model) target space, the slow- roll parameters 
for the spatially fiat Friedmann Universe read (see e.g. [5711115) ): e=-^ {diHf' ~ i|9C/p/L/^ ; rji = 

-if^/iH^i) = 2dlH/H « -e + E;=i3r,C^5jf^/(t^5,[/), where d, d/dip, and |9C/|2 = Y.l=i i.^^Uf. In 
some papers (see e.g. |117H119j ) it was introduced a "cumulative" parameter t] = —Y^^=i'^i'fi/{H\v\^) ~ 
~^ + Yll^ j=i{^i]U){diU){djU) / {U\dU\'^) , where = Yll^=i V'i- One can easily find that for the potential 
l|4.23p parameter rj coincides exactly with parameters r]i and 772. 
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Thus, in these regions it is possible to achieve a period of power-law inflation. 
In spite of a rude character of these estimates, one shall see below that ex- 
ternal space scale factors undergo power-law inflation for trajectories passing 
through these regions. 

Now, let us investigate dynamical behavior of scalar fields and the ex- 
ternal space scale factor in more detail. There are no analytic solutions for 
considered model. So, let us use numerical calculations with applying of a 
Mathematica package proposed in [ jl20] . adjusting it to given models and 
notations (see appendix 3). 

The colored lines on the contour plot of the effective potential in Fig. 14.101 
describe trajectories for scalar fields Lp and with different initial values (the 
colored dots). The time evolution of these scalar fieldqlj is drawn in Fig. 
14. Ill Here, the time t is measured in the Planck times and classical evolution 
starts at t = 1. For given initial conditions, scalar fields approach the global 
minimum of the effective potential along spiral trajectories. 
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Figure 4.11: Dynamical behavior of scalar fields (p (left panel) and 
corresponding initial values denoted by the colored dots in Fig. 14.101 



(right panel) with 



In Figure Iiri2] the evolution of the logarithms of the scale factor a{t) (left 
panel) and the evolution of the Hubble parameter H{t) (right panel) are 
presented and in Fig. 14.131 the evolution of the parameter of acceleration 
q(t). 



"'^^Let us remind that </3 describes fluctuations of the internal space scale factor and reflects the additional 
degree of freedom of the original nonlinear theory. 
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Log[a] H 




Figure 4.12: The number of e-folds (left panel) and the Hubble parameter (right panel) for 
the corresponding trajectories. 




t 



Figure 4.13: The parameter of acceleration (left panel) and its magnification for early times 
(right panel). There are two different form of acceleration with q ^ 1 (De Sitter-like inflation) 
and q ^ 0.75 (power-law inflation with s ~ 4) accordingly. The averaging of q over a 
few periods of oscillations results in g = —0.5 which corresponds to the matter dominated 
decelerating Universe. 
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Because for initial condition the value a{t = 1) = 1 (in the Planck units) is 
used, then log a(t) gives the number of e- folds: log a{t) = N{t). The Figure 
14.121 shows that for considered trajectories one can reach the maximum of 
e-folds of the order of 10. Clearly, 10 e-folds is not sufficient to solve the 
horizon and flatness problems but it can be useful to explain a part of the 
modern CMB data. For example, the Universe inflates by AA^ ^ 4 during the 
period that wavelengths corresponding to the CMB multipoles 2 < / < 100 
cross the Hubble radius jl21Lll22j . However, to have the inflation which is 
long enough for all modes which contribute to the CMB to leave the horizon, 
it is usually supposed that AA^ > 15 fl23] . 

The Figure 14.121 for the evolution of the Hubble parameter (right panel) 
demonstrates that the red, yellow, dark blue and pink lines have a plateau 
H « const. It means that the scale factor a{t) has a stage of the De Sitter 
expansion on these plateaus. Clearly, it happens because these lines reach 
the vicinity of the effective potential saddle point and spend there some time. 

The Fig. 14.131 for the acceleration parameter defined in (IA3.6P confirms 
also the above conclusions. According to Eq. (IA3.8p . g = 1 for the De 
Sitter-like behavior. Indeed, all these 4 lines have stages q ^ 1 for the same 
time intervals when H has a plateau. Additionally, the magnification of this 
picture at early times (the right panel of the Figure 14.131 ) shows that pink, 
green and blue lines have also a period of time when q is approximately 
constant less than one: q « 0.75. In accordance with Eq. (I A3. 81) . it means 
that during this time the scale factor a{t) undergoes the power-law inflation 
a{t) oc with s ~ 4. This result confirms rude estimates made above for the 
trajectories which go through the regions where the effective potential has the 
form (14.231 ). After stages of the inflation, the acceleration parameter starts to 
oscillate. Averaging q over a few periods of oscillations, one obtains q = —0.5. 
Therefore, the scale factor behaves as for the matter dominated Universe: 
a{t) oc Clearly, it corresponds to the times when the trajectories reach 
the vicinity of the effective potential global minimum and start to oscillate 
there. It is worth of noting, that there is no need to plot dynamical behavior 
for the equation of state parameter uj{t) because it is linearly connected with 
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q (see Eq. (IA3.7P ) and its behavior can be easily understood from the pictures 
for q{t). 

As it is shown above for considered quadratic model, the maximal number 
of e-folds is near 10. Can this number be increased? To answer this question, 
let us consider a new model with a higher degree of nonlinearity, i.e. the 
nonlinear quartic model. 



4.3. Quartic model 



In this section let us consider the nonlinear quartic model 

f{R) = R + jR'^- 2Ad . (4.25) 
For this model the scalar field potential (14.161 ) reads [37j : 



5(47)-i/V^-l)'/' + 2Ai^ 



(4.26) 



Here, the scalar curvature R and scalar field are connected as follows: 
e^^ = f' = l + 47i?3 ^ R= [{e^^ - 1)747] . 

Let us look for a solution which has a nonnegative minimum of the effective 
potential Ueff{(p, (j)) (14.131 ) where potential U{4>) is given by Eq. (14.261 ). If 0o 
corresponds to this minimum, then, as it was mentioned above (see also [86]), 
U{(J)o),\ and ff should be positive. To get zero minimum of the effective 
potential, these positive values should satisfy the relation of the form of (14.81 ): 
UicPo) = f? = A/2. 

Additionally, it is important to note that positiveness of U (0o) results in 
positive expression for R{(f)o) > [37]. 





20 


D 





^ -00 =^ [/(0) « ie-^-^ [|(47)-i/3 + 2Ai5] ^ 
+00 and +00 =^ U{4)) ^ |(4^)-i/3e(-s+4A/3)0 _^ ^^^^^^^ 



^"Negative values of and 7 may lead either to negative minima, resulting in asymptotically AdS 
Universe, or to infinitely large negative values of Uef / [37] . In the present consideration let us to avoid both 
of these possibilities. Therefore, one shall consider the case oi Ad,j > 0. See also footnote 
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asymptote -B + 4^/3 = {D - 8)/3^{D - 2){D - 1) < is to be taken 
into account for D < 8. Obviously, the total number of dimensions D = 8 
plays the critical role in quartic nonhnear theories (see [23, EHl [1211 ) 
investigations for D < 8, D = 8 and D > 8 should be performed separately. 
Now, let us consider the case D < 8 (i.e. di = 1, 2, 3). 

It is worth of noting that for considered signs of parameters, the effective 
potential Ueff{<-p,(j)) (14.131 ) acquires negative values when — > +oo (and 
U{(f)) — > 0). For example, if U{(j)Q) = /f = A/2 (the case of zero minimum 
of the effective potential), the effective potential Ueff{(p,(j) ~^ ^) < for 
< e^^'^ < 2 and the lowest negative asymptotic value UefA-^ — > — (16/27)A 
takesplace along the line e^^*^ = 4/3. Therefore, zero minimum of Ue/f is 



21 



local 

As it was mentioned above, extremum positions 0^ of the potential U (0) co- 
incide with extremum positions of Ueff{<-p, (j)) : dU /d(t)\^^ = — > d(f)U(,ff\^. = 
0. The condition of extremum for the potential U (0) reads: 

^.0^i.^-ii±4IiJ + 2A.i^^0. (4.27) 

7(4 -di) 7(4-4) 



For positive 7 and Kr, this equation has two real roots: 



i?o(i) = ^ -W^^^^^-M + VmI , (4.28) 





/ / 2(2 + 4 
2 \^{A-di)k^/M 

where a dimensionless parameter k is introduces as 

k := jAl , (4.30) 

^^It is not difficult to sliow tliat tlie tliin sliell approximation is valid for considered model and a tunnelling 
probability from the zero local minimum to this negative C4// region is negligible. 
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which is positive for positive 7 and A^), and quantities M, u read 

^ ^ _2io/3 (4+^ _ ^ ^ .4 31^ 

3 • 21/3/, (4 - rfi) ' ^^-^^^ 



u = k [-27(4-di)(2 + di)2 

+ Y^272(4 - (ii)2(2 + 4)4 - 4 • 243/c(16 - dlf 



(4.32) 



It can be easily seen that for A; > follows a; < and M > 0. To have real 
w, parameter k should satisfy the following condition 

- 4 • 243(16 - = ■ 

It is not difficult to verify that roots ^0(1,2) are real and positive if < /c < /cq 
and they degenerate for k ^ ko : ^0(1,2) ~^ (A£)/2)a/M. In this hmit the 
minimum and maximum of t/(0) merge into an inflection point. Now, let us 
define which of these roots corresponds to minimum of U (0) and which to 
local maximum. The minimum condition 



defy 



> =^ 7 [{di + 2) - 47i?jj(4 - di)] > (4.34) 



results in the following inequality^ 



7 > : {di + 2)- 47i?g(4 - di) > . (4.35) 



Thus, the root Rq which corresponds to the minimum of U {(/)) should satisfy 
the following condition: 

Numerical analysis shows that ^o(i) satisfies these conditions and corresponds 



to the minimum. For -Ro(2) one obtains that i?o(2) > (^ 4^(4-^^) J cor- 
responds to the local maximum of U{(j)). In what follows one shall use the 



^^As it was already mentioned above, the condition U{(I)q) > leads to the inequality R{4>o) > [37]. 
Taking into account the condition di < 4, it is obvious that inequality (di + 2) + 4|7|_Rg(4 — di) < for 7 < 
cannot be realized. This is an additional argument in favor of positive sign of 7. 
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notations: 

4^min 111 
4^max ^ 111- 



1 + 47^1) , (4.37) 



1 + 47i?, 



3 

0(2) 



(4.38) 



and U{(f)rmn) = Umtn ,U{(f)max) = t^marr- Let US note that (f)rmn , (j^max and the 

ratio Umax/Umin depend on the combination k ( 14.301 ) rather than on 7 and 
Kd taken separately. 

Obviously, because potential ^7(0) has two extrema at and (t)max^ 

the effective potential Ueff{(f,(j)) may have points of extrema only on the 
lines = (J)min and = (l)max where dUeff/d(t)\^^^^^^^^^^ = 0. To find these 
extrema of Uejf-, it is necessary to consider the second extremum condition 
dUeff/dip = on each line separately: 

f -t/„^..-3/fx? + 2Axi = 0, 
^.0^ (4.39) 

[ -/7„,„,-3/2xi + 2Ax2 = 0, 

where xi = exp (^-y/2di/{di + 2) cpi^ > and X2 = exp (^-y/2di/{di + 2) (/?2) 
> 0; and (yi?2 denote positions of extrema on the lines = (j^min and 
4> = (f'max: respectively. These equations have solutions 



Xi(±) = a ± \Ja^ - ^ , 

a > ^=ai] (4.40) 

I On ^m.n.T. 

X2(±) = cx ± \ - p 



« > W/3§^ = «2>«i; (4.41) 

V '-^ min 

where a = A/(3/f) and /3 = Umin/ {?>fl). These equations show that there 
are 5 different possibilities which are listed in the Table [H 

To clarify which of solutions (14.401 ) and (14.411 ) correspond to minima of 
the effective potential (with respect to ^p) one should consider the minimum 
condition 



eff 



>0 =^ /7e,i, + - 4Ax > , (4.42; 
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Table 1: The number of extrema of the effective potential f/g// depending on the relation 
between parameters. 
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where Uextr is either Umin or t/rnaa; and x denotes either xi ov X2- Taking into 
account relations (14.391 ). one obtains 

A 



x'Sft - XA > 



X > 



3/f 



= a. 



(4.43) 



Thus, roots Xi,2(+) define the positions of local minima of the effective poten- 
tial with respect to the variable if and Xi,2(-) correspond to local maxima (in 
the direction of cp). 

Now, let us fix the minimum Xi(+) the point (p = 0. It means that in 
this local minimum the internal space scale factor is stabilized at the present 
day value. In this case 

1 + 13 



Xi(+)|^=o = 1 = a + 



a 



(4.44) 



Obviously, one can do it only io a < 1 ^ G [0,1). For Xi(-) follows: 

Xi(^) = ^• 

Additionally, the local minimum of the effective potential at the point 
((/9 = 0, = (j)min) should play the role of the nonnegative four-dimensional 
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Particular value a — 1 corresponds to the case a = ai = 1 where the only extremum is the inflection 



point with xi{~) = = a = I- Here, A = Umin = 3/f and Ueffi^p = 0, 



= -\ + U„ 
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effective cosmological constant. Tfius, tfie following conditions appear: 

Ae// = Ueff{(p = 0, = (j)rmn) 
= —A + Umin + /l ^ 

=^ -a + (3 + ^>0. (4.45) 

From the latter inequality and equation (14.441) takes place (3 G [|, l)- It 



can be easily seen that P = 1/3 (and, correspondingly, a = 2/3) results in 
Ae// = and one obtains the mentioned above relations: Umin = = -^/2. In 
general, it is possible to demand that A^// coincides with the present day dark 
energy value 10~^^cm~^. However, it leads to very flat local minimum which 
means the decompactification of the internal space [HS]- In what follows, the 
case of zero Ag// is to be mainly considered although all obtained results are 
trivially generalized to Ag// = lO^^^cm^^. 

Summarizing results, in the most interesting case of a > 0^2 the effec- 
tive potential has four extrema: local minimum at ^^/^l^^^^ = 0,0^^^^, lo- 
cal maximum at {(p\x2(-)^ 4'max) and two saddle-points at (^if\x^^_^, (/)min) , and 
('^Ix2(+) '^"^«^) (^^^ ^^S- 14- 151). 



Let us pay particular attention to the case of zero local minimum, when 

^^//('^Ixiw = 0'^—) = ^^^^^ p = 1/3 =^ a = {1 + (3)/2 = 2/3. 
To satisfy the four-extremum condition a > a2, next condition has to be 
satisfied 

U 4 

JY^<L (4.46) 

'-^ min "J 

The fraction Umax /Umin is the function of k and depends parametrically only 
on the internal space dimension di. Inequality (14.461) provides the lower 
bound on k and numerical analysis (see Fig. 14.141) gives k{di = 1) « 
0.000625 ; k{di = 2) « 0.00207 ; k{di = 3) « 0.0035. Therefore, effec- 
tive potentials with zero local minimum will have four extrema if /c G (/c, ko) 
(where /cq is defined by Eq. (14.331) ). The limit k ^ k results in merg- 
ing X2(-) ^ X2{+) and the limit /c — > /cq results in merging Xi(-) ^ X2{-) 
and Xi(+) ^ X2{+)- Such merging results in transformation of correspond- 
ing extrema into inflection points. For example, from Fig. 14.141 follows that 
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Figure 4.14: The form of Umax/Umin as a function of A; G (0, /cq] for di = 1, 2, 3 from left to 
right, respectively. Dashed line corresponds to Umax/Umin = 4/3. 



The typical contour plot of the effective potential with four extrema in the 
case of zero local minimum is drawn in Fig. 14.151 Here, for di = 3 it is taken 
k = 0.004 G {k, ko) which gives a2 « 0.655. Thus, a = 2/3 « 0.666 > a2. 

Let us investigate now a possibility of inflation for considered potential. 
First of all, taking into account the comments in previous section (see Eq. 
( 14.221) ). it is clear that topological inflation in the saddle point Xi(-) well 
as the slow rolling from there in the direction of the local minimum Xi(+) ^ire 
absent. It is not difficult to verified that the generalized power-low inflation 
discussed in the case of the nonlinear quadratic model is also absent here. 
Indeed, from Eqs. (14.131) and (14.261 ) follows that nonlinear potential U (0) can 
play the leading role in the region — > — oo (because U (0) — > for — > +oo). 
In this region Ue/f oc exp {—■\/q<-p) exp (— ^/^0) where q = 2d\/{di + 2) and 
q2 = B"^ = {di + 4)^/[((ii + 2){di + 3)]. For these values of q and q2 the 
slow-roll conditions are not satisfied: e ^ r]i ^ r]2 ^ q/2 + q2/2 > 1. However, 
there are two promising regions where the stage of inflation with subsequent 
stable compactification of the internal space may take place. These are the 
local maximum X2(-) and the saddle X2(+) (see Fig. 14.151 ). Let us estimate 
the slow roll parameters for these regions. 
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Figure 4.15: Contour plot of the effective potential Ueff{^p,(p) (14. 13p with potential f/(0) of 
the form fl4.26p for parameters /3 = 1/3, c^i = 3 and k = 0.004. This plot shows the local 
zero minimum, local maximum and two saddles. The colored lines describe trajectories for 
scalar fields starting at different initial conditions. 

First the local maximum X2(-)- It is obvious that the parameter e is equal 
to zero here. Additionally, from the form of the effective potential (14.131) 
it is clear that the mixed second derivatives are also absent in extremum 
points. Thus, the slow roll parameters r]i and 772, defined in the footnote 
(fTHI). coincide exactly with rjy^ and ?70. In Fig. 14.161 a typical form of these 
parameters is presented as functions of A; G (/c, /co) in the case 13 = 1/3 and 
di = 1,2,3. These plots show that, for considered parameters, the slow roll 
inflation in this region is possible for di = 1,3. 

The vicinity of the saddle point X2(+) is another promising region. Ob- 
viously, if one starts from this point, a test particle will roll mainly along 
direction of (J). That is why it makes sense to draw only Ir/^]. In Fig. 14.171 
typical form of Ir/^l is plotted in the case (3 = 1/3 and di = 1, 2, 3. Left panel 
represents general behavior for the whole range ol k E (A;, /cq) and right panel 
shows detailed behavior in the most interesting region of small k. It shows 
that (ii = 3 is the most promising case in this region. 
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Figure 4.16: Graphs of \ri^\ (left panel) and \ri^\ (right panel) as functions of A; G {k,ko) for 
local maximum X2{-) and parameters (3 = 1/3 and di = 1, 2, 3. 




Figure 4.17: Graphs of Iry^l as functions of k for saddle point X2{+) and parameters ^ = 1/3 
and di = 1,2,3. Left panel demonstrates the whole region of variable k e {k,ko) and right 
panel shows detailed behavior for small k. 
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Now, let us investigate numerically the dynamical behavior of scalar fields 
and the external space scale factor for trajectories which start from the regions 
Xi(-), X2(-) and X2(+)- All numerical calculations perform for /3 = 1/3, di = 3 
and k = 0.004. The colored lines on the contour plot of the effective potential 
in Fig. 14.151 describe trajectories for scalar fields (p and with different 
initial values (the colored dots) in the vicinity of these extrema points. The 
time evolution of these scalar fields is drawn in Fig. 14.181 For given initial 
conditions, scalar fields approach the local minimum Xi(+) of the effective 
potential along the spiral trajectories. 




Figure 4.18: Dynamical behavior of scalar fields (p (left panel) and (j) (riglit panel) witli 
corresponding initial values denoted by the colored dots in Fig. 14.151 

In Figure 14.191 the evolution of the logarithm of the scale factor a{t) (left 
panel) is plotted, which gives directly the number of e-folds and the evolution 
of the Hubble parameter H{t) (right panel) and in Fig. 14.201 the evolution of 
the parameter of acceleration q{t). 

The Figure 14.191 shows that for considered trajectories the maximum of 
e-folds that can be reached is of the order of 22 which is long enough for all 
modes which contribute to the CMB to leave the horizon. 

The Figure 14.191 for the evolution of the Hubble parameter (right panel) 
demonstrates that all lines have plateaus H « const. However, the red, yellow 
and blue lines which pass in the vicinity of the saddle X2(+) have bigger value 
of the Hubble parameter with respect to the dark blue line which starts from 
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Figure 4.20: The parameter of acceleration (left panel) and its magnification for early times 
(right panel). 
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the Xi(-) region. Therefore, the scale factor a{t) has stages of the De Sitter- 
hke expansion corresponding to these plateaus which last approximately from 
100 (dark blue line) up to 800 (red line) Planck times. 

The Fig. 14.201 for the acceleration parameter confirms also the above 
conclusions. All 4 lines have stages q ^ 1 for the same time intervals when H 
has plateaus. After stages of inflation, the acceleration parameter starts to 
oscillate. Averaging q over a few periods of oscillations, one obtains q = —0.5. 
Therefore, the scale factor behaves as for the matter dominated Universe: 
a{t) oc t^/^. Clearly, it corresponds to the times when the trajectories reach 
the vicinity of the effective potential local minimum ^-nd start to oscillate 
there. 

Let us investigate now a possibility of the topological inflation |1251I126] 
if scalar fields cp^cj) stay in the vicinity of the saddle point X2(+)- As it was 
mentioned before, topological inflation in the case of the double-well potential 
takes place if the distance between a minimum and local maximum bigger 
than A(j)cr = 1.65. In this case domain wall is thick enough in comparison 
with the Hubble radius. The critical ratio of the characteristic thickness of 
the wall to the horizon scale in local maximum is r,„ 

ES] and for topological inflation it is necessary to exceed this critical value. 



Therefore, the saddle X2{+) should be checked from the point of these criteria. 

In Fig. 14.211 (left panel) the difference A0 = (J)max — (f'min is drawn for 
the profile ip = V^|x2(+) ^ functions of /c G (/c, /cq) in the case f3 = 1/3 for 
dimensions di = 1,2,3. This picture shows that this difference can exceed 
the critical value if the number of the internal dimensions is c/i = 2 and 
di = 3. Right panel of Fig. 14.211 confirms this conclusion. Here the case 
13 = 1/3, k = 0.004 and rfi = 3 is considered. For chosen values of the 
parameters, A0 = 2.63 which is considerably bigger than the critical value 
1.65 and the ratio of the thickness of the wall to the horizon scale is 1.30 which 
again bigger than the critical value 0.48. Therefore, topological inflation can 
happen for considered model. Moreover, due to quantum fluctuations of 
scalar fields, inflating domain wall will have fractal structure: it will contain 
many other inflating domain walls and each of these domain walls again will 
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contain new inflating domain walls and so on jl25j . Thus, from this point, 
such topological inflation is the eternal one. 



0max-0 
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Figure 4.21: Left panel demonstrates the difference 4>max — (pmin (for the profile (p = V^lx2(+)) 
as a functions of G {k,ko) for parameters /3 = 1/3, and di = 1,2,3 (from left to right 
respectively). Dashed line corresponds to (pmax — 4>min = 1-65. Right panel shows the 
comparison of the potential ^e//(v^lx2{+)' "with a double- well potential for parameters /3 = 
1/3, k = 0.004 and di = 3. 

To conclude this section, let us draw the attention to one interesting fea- 
ture of the given model. From above consideration follows that in the case 
of zero minimum of the effective potential the positions of extrema are fully 
determined by the parameters k and (ii, and for fixed k and di do not depend 
on the choice of A/?. The same takes place for the slow roll parameters. On 
the other hand, if k and di are kept, the hight of the effective potential is 
defined by (see Appendix 2). Therefore, the hight of extrema can be 
changed with the help of A^) but preserve the conditions of inflation for given 
k and di. 

However, the dynamical characteristics of the model (drawn in Figures [4.181 
- 14.201 ) depend on variations of A^) by the self-similar manner. It means that 
the change of hight of the effective potential via transformation Ajj —>■ cAjj 
(c is a constant) with fixed k and di results in rescaling of Figures [4. 181 - 14.201 
in l/\/c times along the time axis. 
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4.4. Bouncing inflation in + model 
4.4.1. The fitting procedure 

Here let us continue the investigation of the quadratic-quartic model started 
in section 2.4 for an arbitrary number of dimensions but in the most par- 
ticular examples = 4 will be used (unless stated otherwise). First of all, let 
us define the relation between the scalar curvature R and the scalaron field 
0. 

Starting by considering a concrete example, for definiteness, let us assume 
7 > 0, a < 0. The pairwise fitting of the various solution branches should 
be performed at points where Q = and different branches of the three- 
solution sector are fitted with each other or to the branches of the one-solution 
sector. The points with Q = correspond to the X— values and X^-, 
Explicitly from (12.611 ) follows 



-max- 



XiR±) = ^«i?± = (4-47) 



and for the concrete configuration 7 > 0, a < 

^max = X{R-) 





1 
















< 
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Xrmn = X {R+) = -a J - — < . (4.48) 



Next, from the defining equation (12.561 ) for the angle that at Q = holds 



so that 



cos(^) = ^ = ^ (4.49) 



Xmax > =^ r>0=^ cos{'d) = 1 

=^ i!) = 27Tm, mGz, 

Xmin < r < =^ cos(i?) = -1 

=^ 7? = TT + 27rn, n G z. (4.50) 
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Figure 4.22: The schematic drawing of the real solution branches and the matching points 
-Pi,2,3,4 This figure shows that points P2,3 (correspondingly, 9 = tc, 2ii) and points Pi^4 (cor- 
respondingly, 6' = 0, Stt) have different nature. So, P2,3 and Pi,4 one shall call branching 
points and monotonic points (in the sense that function R is monotonic in the vicinity of 
these points), respectively. 

Now, the fitting of ttie various solution branclies can be performed as follows 
(see Fig. 14.221) . Let us start with the branch Fi := (i?+ < R < +oo, X^in < 
X < +oo). Moving in on this branch from X ~ +oo one is working in the 
one-solution sector Q > with 



R{ri;Q) = 



r + Q 



1/2 



1/3 



+ 



r-Q 



1/2 



1/3 



(4.51) 



until Q = is hitten at X = X^ax- At this point Pi := (Fi, X = X^ax) £ Fi 
one has to perform the first fitting. Due to r > the following may be chosen 



P(Fi;Q = 0) = 2ri/3 = 2|6| 



1/3 



(4.52) 



so that as simplest parameter choice in (12.581) comes 



Pi = {ri,x = x^ax,Q = o) ^ i9 = o, k = o 

Hence, the parametrization for (Fi, Q < 0) will be given as 

P(Fi,Q < 0) = 2|6|i/3cos(^/3). 



(4.53) 



(4.54) 
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For later convenience, the i!) from the equations ( 12.571 ). ( 12.581 ) is replaced by 
9. The reason will become clear from the subsequent discussion. In should 
be noted that on this Fi— segment one may set "0 = 6. Let us further move 
on Fi until its end at XjrUn^ where again Q = 0. Because there was no other 
point with Q = on this path, the smoothly changing 9 can at this local 
minimum P2 = Fi fl F2 = (X = X^m, R = R+) only take one of the values 
9 = ±7r. For definiteness ^(^2) = tt is chosen. Hence, it holds 

R{P2) = 2|6|i/3cos(7r/3) = \b\^/^ 



= = V-c^/i^l) = R+ (4.55) 

as it should hold. For convenience, let us parameterize the movement on 
the cubic curve by simply further increasing 9. This gives for moving on 
F2 = {R+ > R > R-^X^jiin ^ X < Xj^T^ax) from the local minimum at P2 
to the local maximum at P3 = F2 Pi F3 = [X = X^axiR — R-) ^ further 
increase of ^ by tt up to ^(i^s) = 27r. Accordingly, the complete consistency 
condition looks as 

R{Pi) = 2|6|i/3cos(27r/3) = -|6|i/3 



= -V^=-V-c^/{Ql) = R-- (4.56) 

By further increasing ^ up to 6* = 37r the point P4 = {X = X^m, Q = 0) G F3 
can be reached with 

P(P4) = 2|6|i/3cos(37r/3) = -2\b\^/^ = -2\r\^/\ (4.57) 

Because of r < it can be smoothly fitted it to the one-solution branch 



P(F4, Q)=[r + Q^/'J ' +[r- Q^/^J ' (4.58) 
by setting trivially 

R{P^) = 2{-\r\y/^ = -2\r\^/\ (4.59) 

Summarizing, a very simple and transparent branch fitting picture is achieved, 
where all the movement in the three-solution sector can be parameterized by 
choosing the effective angle as 6* G [0, 37r]. Finally, this picture have to be 



1/3 



1/2 



nl/3 
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fitted in terms of smoottily varying 9 G [0, Stt] with ttie three-solutions Rk 
from (12.581 ). For this purpose let us note that the single value i!) G [0, tt] in 
( 12.581 ) is a projection of the smoothly varying 6 G [0, Stt]. Fixing an arbitrary 
'd one easily finds the following correspondences 

^(Fi, ^) = e{T2, 1?) = 27r - ^(Fs, d) = 2Ti + d (4.60) 

and hence 

RWiJ)] = 2|6|i/3cosf|') =%.=o) = ^3, 



R[e{r2,i^)] = 2|6|i/^cos 



3 




R[eir,,^)] = 2|6|i/3cos, ^ 

= %=l) = 7^l. (4.61) 

Analogically, the rules for fitting procedure can be obtained in the case 7 < 
0, Of > 0. So, all the fitting mechanism is clear now and can be used in further 
considerations. 

4.4.2. Dynamics of the Universe and scalaron in R^ + R^ model 

To study the dynamics of the Universe in given model, let us assume that 
the four-dimensional metric g in (12.61 ) is spatially flat Friedmann- Robertson- 
Walker one: 

g = -dt^dt + a^{t)dx^dx. (4.62) 

Thus, scalar curvatures R and R and scalaron are functions of time. There- 
fore, Eq. (12.351 ) for homogeneous field reads 

dU 

+ 3//0 + — = 0, (4.63) 
a(p 

where the Hubble parameter H = d/a and the dotes denote the differentiation 
with respect to time t. Potential U is defined by Eq. (12.51 ). Because U 
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Figure 4.23: The form of the potential fl3.8p as a multi- valued function of X = e^'^ — 1 in 
the case D = 4, = 0.1, 7 = 1 and a = —1. Points -Pi, 2,3,4 are defined in Fig. 14.221 



depends on R which is a multi-valued function of (f) (or, equivalently, o 



X), the potential U is also a multi- valued function of X (see Fig. 14.231) 
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However, previous analysis shows that this problem can be avoided by making 
X and R single- valued functions of a new field 6 (one would remind that the 
particular case of a7 < is considered when a < 0, 7 > 0): 



x{e) = { 



^(1/7) (2|«|/3)' cosh(^) ,^<0 
^(1/7) (2|a|/3)^ cos(^) , < ^ < 37r 
^(1/7) (2|a|/3)^ cosh(^ - 37r) , ^ > 37r 



(4.64) 



and 



2v^|a|/(67) cosh(^/3) , ^ < 

2^\a\/{6-f) cos{e/3) , < ^ < 37r 
-2v^|a|/(67) cosh[(^/3) - tt] , 6 > Stt 



(4.65) 



The function X = X(6) is schematically given in Fig. 14.241 It is necessary 
to keep in mind that the case under consideration is /' > — > X > — 1. If 
condition X^in > — 1 is imposed (in opposite case the graphic X{6) will be 
cut into two disconnected parts) then the parameters a and 7 should satisfy 



^^In spite of the divergency of dR/dX in the branching points -Pj.a, the derivatives dU/dX are finite in 
these points. Moreover, R and X have the same values in branching points for different branches. Therefore, 
the branches arrive at the branching points with the same vahies of dU/dX and Fig. 14.231 clearlv shows it. 
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X = X„,„ cosh(^) X = X„„ cos(0) X = X,„.„ cosh(6l - 3;r) 

0<O Q<0<1,7T S>^K 



Figure 4.24: The schematic drawing of Eq. fl4.64p in the case X^j^ > —1. Here, X„ 
^(1/7) (2|«|/3)% = -y/(l/7) (2|a|/3f and B^^. is defined by Eq. (USB. 



the inequality: 



X,„i„ = -^(1/7) (2|a|/3f > -1 



'^-,1/3 



(4.66) 



The maximal value of Q (which is greater than Stt in the case Xj^i^ > ~1) is 
defined from the transcendental equation 



where c := 



1 
2 

= cosh[(^™aa;/3) - tt] , 
\ -1 



1/3" 



(4.67) 



1/7) (2|a|/3)^ 



The hmit X 



1 corresponds to the 



limit 6 — > 6 max- With the help of Eq. (14.641 ) and formula 

1 dX 



de A{X^l)de 
the following useful expressions can be also obtained: 

dX 



(4.68) 



de 



:0,7r,27r,37r 



de 







(4.69) 



6'=0,7r,27r,37r 



It can be easily verified that field e satisfies the equation: 



e + sm + r{e)e'^ + G{e) 



dU 







(4.70) 
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Here, G{6) = G^^ = (Gii) ^ = {dcjy/dO) ^ is the one dimensional metric on 
the moduh space with the corresponding Christoffel symbol r(^) = Pj^ = 

(1/2)^11(^11), = {£(j)/de^)/{d(j)/de). 

4.4.3. Properties of the potential U{9) 

As it was mentioned above, in current case the potential (12.51) as a function 
of ^ is a single- valued one. Now, let us investigate analytically some general 
properties of U{6). In this section, D is an arbitrary number of dimensions 
and signs of a and 7 are not fixed if it is not specified particularly. 

First, let us consider the extrema of the potential U{6). To find the ex- 

tremum points the equation have to be solved 

dU dUd(j) dU 1 dX ^ 

— = - = = 0. (4.71) 

de d(i)de d(i)A{x + i)de ^ ' 

Therefore, the extrema correspond either to the solutions of the equation 
( lOm dU/d(t) = for finite dX/dO (X > -1) or to the solutions of the 
equation dX/d9 = (X > — 1) for finite dU/d(j). The form of the potential U 
(see Eq. (12.51) ) shows that this potential and its derivative dU/d(j) is finite for 
X > —1. Thus, as it follows from Eq. (14.691) . the potential U{9) has extrema 
at the matching points ^ = 0, tt, 27r, Stt. Additional extremum points are real 
solutions of the equation (12.691 ). For the given model (12.431 ) this equation 
reads: 

The form of this equation shows that there are two particular cases: D = 8 
and D = 4. The D = 8 case was considered previously. Let us consider now 
the case D = 

R'--R+^ = 0. (4.73) 
27 7 

It is worth of noting that parameter a disappeared from this equation. Thus 
a has no an influence on a number of additional extremum points. To solve 
this quartic equation, one should consider an auxiliary cubic equation 

o 8A4 1 , , 

-u r = 0. 4.74 

7 472 ^ ^ 
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The analysis of this equation can be performed in similar manner as it was 
done it for the cubic equation (12.451 ). Let us introduce the notations: 



q — f — ^ ^ 



3 7 



'7' 
1 f 1 



'8A4 



T 



82 



(4.75) 



It make sense to consider two separate cases. 
1. sign 7 = — sign A4 =^ Q > . 

In this case only one real solution of Eq. (14.741) is needed: 



Ui 



1/2 



1/3 



+ 



f-Q 



1/2 



1/3 



> 



(4.76) 



Then, solutions of the quartic (14.731) are the real roots of two quadratic equa- 
tions 



E? ± y/uiR + ^ (uizt e\l + 3^ ) = , e = sign 7 



(4.77) 



Simple analysis shows that for any sign of 7 two real solutions can be obtained: 

7 < ^ = -lu'f ± + 1(^2 + 3^)1/2, 



(4.78) 



7 



> ^ = \u^f ± ^-iwi + i(w2^3g)V2 



2. sign 7 = sign A4 ^ Q ^ . 

It is not difficult to show that in this case the real solutions of the form 
of (14.781) (where the evident substitution is done u{ + ?>q ^ Ui — ?>\q\) takes 
place if 



g>o 



1 1/3 



< 



32IA4 



Now, let us investigate zeros of the potential U. For f'y^O 
the condition of zeros of the potential (12.51) is: 

Rf - f = + aR^ + 2Az? = . 

Therefore, zeros are defined by equation: 

2 



(4.79) 
(4.80) 



a 



R^ = -—± 
07 



a 
67, 



2A 



D 



37 



1/2 



(4.81) 
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Obviously, the necessary conditions for zeros are: 

7 > => Ad< a7(247) , 



(4.82) 



7 < 



> -«V(24|7|) 



Additionally, the positiveness of the r.h.s. of the equation (14.811) should be 
checked. 

Let us consider now asymptotical behavior of the potential U{0). Here, 
the limits 6 — > 6max and 6 — )• — oo. In the former case it appears: 



max 



U{6) -sign [f [9 max)) X oo 



(4.83) 



In the latter case: 



9 



-oo 



U{e) - exp 9\ ^ 



+00 



D > 



const > ,D = 
+0, 2<D< 



(4.84) 



where Eqs. (14.641 ) and (14.651) have been used. This result coincides with 
conclusions of appendix 1 in 

To illustrate the described above properties, let us draw the potential 
U{9) in Fig. 14.251 for the following parameters: D = 4, A4 = 0.1, 7 = 1 and 
a = —1. These parameters contradict to the inequahties (14.791) and (14.821) . 
Therefore, 9 = 0, tt, 27r, Stt are the only extremum points of the potential 
U{9) and zeros are absent. These parameters are also satisfy the condition 
(14.661 ). The absence of zeros means that all minima of the potential U{9) are 
positive. 

For the subsequent investigations, it is useful also to consider an effective 
force and mass squared of the field 9. As it follows from Eq. (14.701 ). the 
effective force: 



F = -G{9) 



dU 

1b 



(4.85) 
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Figure 4.25: The form of the potential f l2.5p as a function of ^ in the case D = 4, A4 = 0.1, 
7 = 1 and a = —1. For these values of the parameters, all extrema correspond to the 
matching points = 0, tt, 27r, Stt. In the branching points ^ = tt, 27r the potential has local 
maximum and local non-zero minimum, respectively, and the monotonic points 6' = 0, Svr are 
the inflection ones. Potential tends asymptotically to +00 when Q goes to Qraax and to zero 
when d — 7- —00. 



Varying Eq. ( 14.701 ) with respect to field Q one obtains dynamical equation 
for small fluctuation bQ where mass squared reads: 



m\ = G{e) 



d?U dG{e) dU 



+ 



(4.86) 



' de de 

In Fig. I4.26l the effective force and the mass squared are presented as functions 
of 9 for the potential drawn in Fig. 14.251 These figures indicate that field 
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Figure 4.26: The effective force fl4.85p (left panel) and the mass squared (14.861) (right panel) 
for the potential U{9) drawn in Fig. 14.251 These pictures clearly show singular behavior of 
F and in the matching points 6' = 0, tt, 27r, Stt. 



9 may have very nontrivial behavior. This non-triviality follows from two 
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reasons. First, the field 6 has non-canonical kinetic term which result in 
appearing of non-flat moduli space metric G{6) and derivative of G{6) in 
Eq. (14.701 ). Second, the function G{9) has singular behavior at the matching 
points ^ = 0, TT, 27r, Stt. Thus, the given intuition does not work when one 
wants to predict dynamical behavior of fields with equations of the form of 
(14.701 ) with potential drawn in Fig. 14.251 especially when field approaches 
the matching points. It is necessary to solve equations analytically or to 
investigate them numerically. Such analysis for the given model is performed 
in the next section where the concentration is focused on the case, where all 
extrema correspond to the matching points ^ = 0, tt, 27r, Stt. 

4.4.4. Dynamical behavior of the Universe and field 9 

Let us investigate dynamical behavior of scalar field 9 and the scale factor 
a in more detail. There are no analytic solutions for considered model. So, 
numerical calculations are used. To do it, a Mathematica package proposed 
1120] is applied and adjusted to the given models and notations in appendix 



m 



3 of the paper According to these notations, all dimensional quantities in 
the graphics are given in the Planck units. Additionally, one should remember 
that metric on the moduli space is not flat and defined in Eq. (14.701 ). For 
example, the canonical momenta and the kinetic energy read: 



Pe = -Gii9 = -[^\ 9, (4.87) 

t^i Kl \d9 J 

where 87rG = ti\ and G is four-dimensional Newton constant. To understand 
the dynamics of the Universe, one shall also draw the Hubble parameter: 

3 f-^ =?>H'' = -Gii9^ + Ui9) (4.88) 
\a 2 



and the acceleration parameter: 
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Fig. 14.261 shows that the effective force changes its sign and the mass 
squared preserves the sign when 6 crosses the branching points tt, 27r and vise 
verse, the effective force preserves the sign and the mass squared changes the 
sign when 9 crosses the monotonic points 0, Stt. Therefore, it make sense to 
consider these cases separately. 



4.4.5. Branching points ^ = tt, 27r 



First, the dynamical behavior of the Universe and scalaron is considered 
in the vicinity of the branching point 6 = 2ti which is the local minimum 
of the potential in Fig. 14.251 The time evolution of scalaron field 9 and its 
kinetic energy Ekin are drawn in Fig. 14.271 Here and in all pictures below, 
the same parameters are used as in Fig. 14.251 The time t is measured in 
the Planck times and classical evolution starts at t = 1. For the initial 
value of 9 it is taken that 9initiai = 3.5. In Fig. 14.281 the evolution of the 




Figure 4.27: Dynamical behavior of scalar field d{t) (left panel) and its kinetic energy Ekin{t) 
(right panel) in the vicinity of the branching point 9 = 2tt. 



logarithms of the scale factor a{t) (left panel) and the evolution of the Hubble 
parameter H{t) (right panel) is plotted, and in Fig. 14. 291 the evolution of the 
parameter of acceleration q{t) (left panel) and the equation of state parameter 
u{t) = [2q{t) + l]/3 (right panel). 

Fig. 14.271 demonstrates that scalar field 9 bounces an infinite number 
of times with decreasing amplitude in the vicinity of the branching point 
9 = 27T. 9 cannot cross this point. From Figs. 14.281 and 14.291 it can be 
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Figure 4.28: The time evolution of the logarithms of the scale factor a{t) (left panel) and 
the Hubble parameter H{t) (right panel) for the trajectory drawn in Fig. 14.271 A slightly 
visible oscillations of In (a) (caused by bounces) can be seen by magnification of this picture 
(middle panel). 




Figure 4.29: The parameter of acceleration q{t) (left panel) and the equation of state pa- 
rameter u{t) (right panel) for the scale factor in Fig. 14.281 
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seen that the Universe asymptotically approaches the de Sitter stage: H — ^ 
const, q — > +1, and u — > —1. Such accelerating behavior can be called 
bouncing inflation. 

Concerning the dynamical behavior in the vicinity of the branching point 
9 = TT^ the analysis (similar performed above) shows that the scalaron field 
9 cannot cross this local maximum regardless of the magnitude of initial 
velocity in the direction of ^ = tt. It bounces back from this point. 



4.4.6. Monotonic points ^ = 0, Stt 



Now, let us investigate the dynamical behavior of the model in the vicinity 
of the monotonic points ^ = 0, Stt which are the points of inflection of the 
potential in Fig. 14.251 Figs. 14.251 and 14.261 show that for both of these points 
the model has the similar dynamical behavior. Therefore, for definiteness, 
let us consider the point 9 = Stt. To investigate numerically the crossing of 
the monotonic point Stt, it is necessary to take very small value of a step 
At. It can be achieved if very large value of the maximum number of steps is 
chosen. Thus, for the given value of the maximum number of steps, the closer 
to Stt the initial value 9initiai is taken the smaller step At is obtained. For 
current calculation let us choose 9initiai = 9.513. Fig. 14.301 demonstrates that 
scalar field 9 slowly crosses the monotonic point 37r with nearly zero kinetic 
energjo- Then, just after the crossing, the kinetic energy has its maximum 
value and starts to decrees gradually when 9 moves to the direction 27r. 
Figs. |4,31J and [432] demonstrate the behavior of the Universe before and after 
crossing 37r. The vicinity of the branching point 27r is not shown because 
when 9 approaches 27r the Universe has the bouncing inflation described 
above. Hence, there are 3 phases sequentially: the short de Sitter-like stage 
during slow rolling in the vicinity of the inflection point before crossing, then 
decelerating expansion just after the crossing with gradual transition to the 
accelerating stage again when 9 approaches the branching point 27r. Clearly 



^^The derivative dO/dt goes to — oo when 6 ^ in (with different speed on different sides of Stt) but 
: at Stt and kinetic energy is finite (see Eq. (14.871) ). 
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Figure 4.30: Dynamical behavior of scalar field 6(t) (left panel) and its time derivative 
d9/dt (middle panel) and kinetic energy Ekmit) (right panel) for the case of crossing of the 
inflection point 9 = 3n. 
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Figure 4.31: The time evolution of the logarithms of the scale factor a{t) (left panel) and 
the Hubble parameter H{t) (right panel) for the trajectory drawn in Fig. 14.301 




Figure 4.32: The parameter of acceleration q{t) (left panel) and the equation of state pa- 
rameter u(t) (right panel) for the scale factor in Fig. 14.311 
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Figure 4.33: The time evolution of the logarithms of the scale factor a{t) (left panel) and the 
Hubble parameter H{t) (right panel) for the trajectory drawn in Fig. I4.27l in the Brans-Dicke 
frame. 

that for another monotonic point ^ = the similar crossing behavior (without 
the bouncing stage when 6 — > — oo) can be reached. Therefore, the monotonic 
points ^ = and 9 = Sir are penetrable for the scalaron field 9. 

In appendix 4, it is shown that bouncing inflation in the vicinity of the 
branching point takes place also in the Brans-Dicke frame. 

Summary 

A possibility of inflation in multidimensional cosmological models has been 
investigated. The main attention was paid to nonlinear (in scalar curvature) 
models with quadratic and quartic lagrangians. These models con- 
tain two scalar fields. One of them corresponds to the scale factor of the 
internal space and another one is related with the nonlinearity of the original 
models. The effective four-dimensional potentials in these models are fully 
determined by the geometry and matter content of the models. The geometry 
is defined by the direct product of the Ricci-flat external and internal spaces. 
As a matter source, a monopole form field, D-dimensional bare cosmological 
constant and tensions of branes located in fixed points have been included 
into consideration. The exact form of the effective potentials depends on the 
relation between parameters of the models and can take rather complicated 
view with a number of extrema points. 

First of all, it was found that a range of parameters which insures the 
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existence of zero minima of the effective potentials. These minima provide 
sufficient condition to stabilize the internal space and, consequently, to avoid 
the problem of the fundamental constant variation. Zero minima correspond 
to the zero effective four-dimensional cosmological constant. In general, one 
can also consider positive effective cosmological constant which corresponds 
to the observable now dark energy. However, it usually requires extreme fine 
tuning of parameters of models. 

Then, for corresponding effective potentials, the possibility of the external 
space inflation was investigated. It have been shown that for some initial 
conditions in the quadratic and quartic models can be achieve up to 10 and 
22 e-folds, respectively. An additionally bonus of the considered model is 
that model can provide conditions for the eternal topological inflation. 

Obviously, 10 and 22 e-folds are not sufficient to solve the homogeneity 
and isotropy problem but big enough to explain the recent CMB data. To 
have the inflation which is long enough for modes which contribute to the 
CMB, it is usually supposed that AN > 15 [I23j. Moreover, 22 e-folds is 
rather big number to encourage the following investigations of the nonlinear 
multidimensional models to flnd theories where this number will approach 
50-60. It can be seen that the increase of the nonlinearity (from quadratic to 
quartic one) results in the increase of AA^ in more that two times. So, there is 
a hope that more complicated nonlinear models can provide necessary 50-60 e- 
folds. Besides, this number is reduced in models where long matter dominated 
(MD) stage followed inflation can subsequently decay into radiation jl27Lll28] . 
Precisely this scenario takes place for given models. It also have been shown 
for quadratic and quartic nonlinear models, that MD stage with the external 
scale factor a ~ t'^^^ takes place after the stage of inflation. It happens 
when scalar flelds start to oscillate near the position of zero minimum of the 
effective potential. 

For the investigation of the dynamical behavior of the scalaron fleld 
and the Universe in nonlinear model with curvature-squared and curvature- 
quartic correction terms: , the parameters a and 7 are chosen in such a way 
that the scalaron potential ^7(0) is a multi- valued function consisting of a 
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number of branches. These branches are fitted with each other either in the 
branching points (points P2,3 in Fig. 14.231 ) or in the monotonic points (points 
Pi^4 in Fig. 14.231 ). The potential U was optimized so that a way that it 
becomes the one- valued function of a new field variable 6 = 6 {(J)) (see Fig. 
14.251 ). This has enabled to consider the dynamical behavior of the system in 
the vicinity of the branching and monotonic points in [D = 4)-dimensional 
space-time. This investigation shows that the monotonic points are pene- 
trable for scalaron field (see Figs. I4.30M.32I ) while in the vicinity of the 
branching points scalaron has the bouncing behavior and cannot cross these 
points. Moreover, there are branching points where scalaron bounces an infi- 
nite number of times with decreasing amplitude and the Universe asymptot- 
ically approaches the de Sitter stage (see Figs. I4.27ti4.29l ). Such accelerating 
behavior is called bouncing inflation. It should be noted that for this type 
of inflation there is no need for original potential U (0) to have a minimum 
or to check the slow-roll conditions. A necessary condition is the existence 
of the branching points. This is a new type of inflation. It is shown that 
this inflation takes place both in the Einstein and Brans-Dicke frames. It is 
found that this type of inflation for the model with the curvature-squared 
and curvature-quartic correction terms which play an important role during 
the early stages of the Universe evolution. However, the branching points 
take also place in models with R^^-type correction terms jl29j . These terms 
play an important role at late times of the evolution of the Universe. 

There is no need for fine tuning of the initial conditions to get the bouncing 
inflation. In Figs. (I4.27I )-( H.29I ). for definiteness the initial conditions are 
chosen as ^ = 3.5 and Ekin = 0. However, the calculations show that these 
figures do not qualitatively change if one takes arbitrary 6 G (tt, 27t) and 
non-zero E^in. Fig. [4,261 indicates that the minimum at ^ = 27r is stable with 
respect to tunneling through the barrier at this point. The situation is similar 
to the quantum mechanical problem with infinitely high barrier. It has been 
stressed already that the form of the potential Fig. 14.251 is not sufficient to 
predict the dynamical behavior of 6. This field has very non-trivial behavior 
because of non-canonical kinetic term and singular (at the matching points) 
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non-flat moduli space metric G{6). Therefore, it is impossible to "jump" 
quantum mechanically from one branch to another. It is impossible to apply 
to given dynamical system the standard tunneling approach (e.g., in jl30j ). 
This problem needs a separate investigation. It is worth of noting that the 
Universe with a bounce preceding the inflationary period was considered in 
1131] where it was shown that due to a bounce the spectrum of primordial 



perturbations has the characteristic features. It indicates that the similar 
effect can take place in this model. 
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CONCLUSION 

Several nonlinear multidimensional gravitational models are studied in this 
Thesis from the point of freezing stabilization and compactification of ad- 
ditional dimensions; and possibility of inflation. The main results may be 
constituted as follows. 

Two separate branches of the four-dimensional spaces are produced for 
multidimensional nonlinear models with the scalar curvature of type 
after dimensional reduction. For these type of models, simultaneous com- 
pactification and late-time acceleration is unreachable, due to branching of 
potential. 

The pure gravitational model with curvature-quadratic and curvature- 
quartic correction terms allows too the stable compactification of the internal 
space for certain interval of the parameters (regions of stability). Assuming 
the scale factor of the stabilized internal space to be in order of the Fermi 
length: 6(o)i ^ Lp ^ lO^^^cm, one arrive that a ~ Ljp and for the effective 
cosmological constant and masses folds: — Ae// ~ ~ ~ ITeV^. 

Adding forms to the model of ^^^gives possibility to obtain positive mini- 
mum of the effective potential. In this model the freezing stabilization of the 
internal spaces is provided. This allows to avoid the problem of the funda- 
mental constant variation in multidimensional models [251196]. At the same 
time, the stage of the cosmic acceleration is achieved too. However, fine tun- 
ing is required to get parameters of the present-day accelerating expansion. 

Additionally the existence of domain walls was found for this model, 
which separates regions with different vacua in the Universe without pro- 
viding of inflation, because the effective potential is not flat enough around 
the saddle point. In case of parameter z < 2 the minimums of the poten- 
tial is metastable, i.e. the quantum tunneling is possible both in and in 
if directions. For z > 2, tunneling in ip direction is still possible due to 
Ueff{(f,(J)o) « e^'^U{(J)o) — > for (yi) — > — oo which is less than any positive 
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Ag//. This could result in the materialization of bubbles of the new phase in 
the metastable one (see e.g. [27]). Hence, late-time acceleration is possible 
only if characteristic lifetime of the metastable stage is greater than the age 
of the Universe. For this model the global negative minimum is also possible. 
This minimum is stable both in classical and quantum limits. However, the 
acceleration is absent because of its negativity. 

Another feature of this model consists in multi- valued form of the effective 
potential, which was two branches for for each choice of parameter /i. This 
gives possibility for investigation of transitions from one branch to another 
one by analogy with catastrophe theory or similar to the phase transitions 
in statistical theory. In this particular model the point = corresponds to 
the singularity R, R ^ ±oo. Indicating that the analog of the second order 
smooth phase transition through the point = is impossible in this model. 
The analog of the first order transition via quantum jumps is possible too. 

The inflation has been investigated in the models with nonlinearities of 
type R^ and R'^ with a monopole form field, D-dimensional bare cosmological 
constant and tensions of branes located in fixed points as source. As result, 
the equivalent linear model has two scalar fields: one corresponds to the scale 
factor of the internal space and another the the nonlinearity of the original 
models. Up to 10 and 22 e-folds can be achieved for R'^ and respectively, 
although it is not sufficient to solve the problems of homogeneity and isotropy 
in recent CMB data, where 50-60 e-folds are needed. This number is reduced 
in case if inflation, followed by long matter dominated stage, can subsequently 
decay into radiation |127l[T28] . This scenario takes place for given models, 
as have been shown for R^ and models: matter dominated stage with the 
external scale factor a ~ 

^2/3 

occurs after the stage of inflation. It happens 
when scalar fields start to oscillate near the position of zero minimum of the 
effective potential. Another significant problem for these models consists in 
the spectral index. For example, in the case of model, the spectral index: 
Ug ~ 1 + 2?7|^2(+) ~ 0-61 which is less than observable now ^ 1. As possible 
solution of this problem, more complicated function f{R) can be used. For 
example, in |t5B] it is shown that quadratic and quartic nonlinearities lead to 
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increase of Ug. 

As next step, the dynamical behavior of the scalaron field (f) and the Uni- 
verse in quadratic and quartic nonlinear model was investigated. In this 
case, the scalaron potential U{(j)) can be a multi-valued function consisting 
of a number of branches, being connected in the branching points or in the 
monotonic points. Further reparametrization of the potential U to the one- 
valued function of a new field variable 6 = 9[(j))^ enables consideration of the 
dynamical behavior of the system in the vicinity of the branching and mono- 
tonic points. It appears that the monotonic points are penetrable for scalaron 
field while in the vicinity of the branching points scalaron has the bouncing 
behavior and cannot cross these points. Also branching points where found, 
where scalaron bounces an infinite number of times with decreasing ampli- 
tude, approaching asymptotically de Sitter stage. This king of accelerating 
behavior is called bouncing inflation. For this type of inflation there is no 
need for original potential U (0) to have a minimum or to check the slow-roll 
conditions, i.e. no need for fine tuning. The existence of the branching points 
is only condition imposed. The minimum at ^ = 27r is stable for tunneling 
through the barrier at this point. 

This represents a new type of inflation, which takes place both in the 
Einstein and Brans-Dicke frames and play an important role during the early 
stages of the Universe evolution. In the model of nonlinearity [ ll29j the 
branching points also take place, playing an important role at late times of the 
evolution of the Universe. Therefore, bouncing inflation may be responsible 
for the late-time accelerating expansion of the Universe. 
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APPENDIXES 



Appendix 1. The critical dimensions in f{R) theories 

The existence of a critical dimension (in this case D = 8) is a rather 
general feature of gravitational theories with polynomial scalar curvature 
terms (see, e.g., Refs. |i21ll32l[T^ ). Following the paper [37], it can be easily 
demonstrated for a model with curvature nonlinearity of the type 

TV 

fiR) = J2akR' (Al.l) 

k=0 

for which the ansatz 

N 

e^^^ = f = Y^kakR^-^ (A1.2) 

fc=0 

leads, similar like (12.51 ). to a potential 

1 ^ 
/7(0) = - (/)-^/(^-2) ^(/c - l)akRK (A1.3) 



The condition of extremum (12.301 ) for this potential reads 



N 

Df - 2Rf = ^ (L» - 2k) akR^ = . (A1.4) 

fc=0 

Thus, at the critical dimension D = 2A^ the degree of this equation is reduced 
from to A^ — 1. In this case the search of extrema is considerably simplified. 

In the limit — > +oo the curvature will behave like R ~ ce^'^ where h and 
c can be defined from the dominant term in (IA1.2P : 

e^^ « NaNR""-^ - Aa^c^'^e^^^i)'^'^. (A1.5) 
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Here the requirement /' > allows for the following sign combinations of the 
coefhcients and the curvature asymptotics R{(t) — > oo): 



N = 21 : sign [a^] = sign — > oo^ 

iV = 2/ + l: aiv>0, sign ^ oo)] = ±1. (A1.6) 

The other combinations, N = 21 : signfa^r] = — sign [^(0 — > oo)], = 
2/ + 1 : Oat < 0, sign [R{(j) — > oo)] = ±1, would necessarily correspond to 
the f < sector, so that the complete consideration should be performed in 
terms of the extended conformal transformation technique of Ref. jSD] • Such 
a consideration is out of the scope currently and the focus is restricted to the 
cases (I A 1.61) . The coefficients h and c are then easily derived as h = A/ (N—l) 



and c = sign {un) lAa^vl Plugging this into (IA1.3P one obtains 

^ +oo) « sign I AaA^r^ e^^^'^e^^'^ (A1.7) 
and that the exponent 

D-2N ^ , 

-A (A1.8) 



{D-2){N-r 
changes its sign at the critical dimension D = 2N: 

Ioo for L» > 2 A , 
1 for L» = 2 A , 
for L> < 2 A . 

(A1.9) 

This critical dimension D = 2N is independent of the concrete coefficient 
ttN and is only defined by the degree deg^(/) of the scalar curvature polyno- 
mial /. From the asymptotics ( lAl.Op follows that in the high curvature limit 
(J) —>■ +00, within given oversimplified classical framework, the potential U{(f)) 
of the considered toy-model shows asymptotical freedom for subcritical di- 
mensions D < 2 A, a stable behavior for a^v > 0, D > 2 A and a catastrophic 
instability for a^v < 0, D > 2 A. It should be noted that this general behavior 
suggests a way how to cure a pathological (catastrophic) behavior of poly- 
nomial ^^1— nonlinear theories in a fixed dimension D > 2Ni. By including 
higher order corrections up to order A2 > D/2 the theory gets shifted into 
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the non-pathological sector with asymptotical freedom. More generally, one 
is even led to conjecture that the partially pathological behavior of models in 
supercritical dimensions could be an artifact of a polynomial truncation of an 
(presently unknown) underlying non-polynomial f{R) structure at high cur- 
vatures — which probably will find its resolution in a strong coupling regime 
of M— theory or in loop quantum gravity. 



Appendix 2. Self-similarity condition 



Due to the zero minimum conditions U{(j)rnin) = fi = A/2, the effective 
potential (14.131 ) can be written in the form: 



2di 



Ueff{(f, 0) = U{(j)rmn)e 

(A2.1 



X 



[/(0) ^ ridV _ rid 



U{(j)jn, 

Exact expressions for U (0) (14.201 ) and (14.261 ) indicate that the ratio 

-^ = F(0,/c,rfi) (A2.2) 

depends only on 0, /c and di. Dimensionless parameter k = ^Aj:) for the 
quadratic model and k = 7A|, for the quartic model. In Eq. (A2.2p one have 



to take into account that ^^mi is a function of k and di: (J)min = 0rnm(^, o^i)- 
Then, t/((/)^,„) defined in Eqs. ( Km and ( Km reads: 



U{(t)rmn) = Ai?F(0min(A:, C^i) , /c, c/i) . (A2.3) 

Therefore, parameters k and di determine fully the shape of the effective po- 
tential, and parameter Aj:) serves for conformal transformation of this shape. 
This conclusion is confirmed with the fact that all extrema depend only on k 
and di. Thus, figures I^TTUl and 14.151 for contour plots are defined by k and di 
and will not change with Ajj. From the definition of the slow roll parameters 
it is clear that they also do not depend on the hight of potentials and in the 
given model depend only on k and di (see figures 14.161 and 14.471 ). Similar 
dependence takes place for difference A0 = ^^aa; ~ 0mm drawn in Fig. 14.211 
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Thus the conclusions concerning the slow roll and topological inflations are 
fully determined by the choice of k and di and do not depend on the hight 
of the effective potential, in other words, on Ajj. So, for fixed k and di pa- 
rameter Ad can be arbitrary. For example, one can take Ajj in such a way 
that the hight of the saddle point X2(+) will correspond to the restriction on 
the slow roll inflation potential (see e.g. fl34j ) Ueff < 2.2 x lO^^^Mp^, or in 
the notations t/g// < 5.5 x lO^^'^Mj,;. 

Above, the figure is indicated which (for given k and di) do not depend on 
the hight of the effective potential (on Ad). What will happen with dynamical 
characteristics drawn in figures STTSl 14A9| and [4,201 ( and analogous ones for the 
quadratic model) if keeping fixed k and di, will change A^)? In other words, 
the positions of the extrema points (in ((/?, 0)-plane)are kept, but the hight 
of extrema experience change. One can easily answer this question using the 
self-similarity condition of the Friedmann equations. Let the potential U in 
Eqs. (IA3.2P and (I A3. 31} be transformed conformally: U ^ cU where c is a 
constant. Next, let us introduce a new time variable r := -y/ct. Then, from 
Eqs. ( IA3.2p -( lA3.5p follows that the Friedmann equations have the same form 
as for the model with potential U where time t is replaced by time r. This 
condition one may call the self-similarity. Thus, if in the model one changes 
the parameter A^ : A/j— ^cA/j, it results (for fixed k and di) in rescaling 
of all dynamical graphics (e.g. Figures 14.181 - 14.201 ) along the time axis in 
l/^/c times (the decrease of A/j leads to the stretch of these figures along 
the time axis and vice versa the increase of Ajj results in the shrink of these 
graphics). Numerical calculations confirm this conclusion. The property of 
the conformal transformation of the shape of Ueff with change of A/j for fixed 
k and di can be also called as the self-similarity. 



Appendix 3. Multi-component scalar field model 



Let us consider n scalar fields minimally coupled to gravity in four dimen- 
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sions. The effective action of this model reads 

- Gi0'<'>'"'d^y^%^-2U{v\v^...)J (A3.1) 

where the kinetic term is taken in the canonical form: Gij = diag(l, 1, . . .) 
(flat a model). Such multi-component scalar fields originate naturally in 
multidimensional cosmological models (with linear or nonlinear gravitational 
actions) |571lll0l[TT5] . Here, the the usual conventions c = h = 1 is used, i.e. 
Lpi = tpi = 1/Mpi and SttG = Sir/Mpi. Here, scalar fields are dimensionless 
and potential U has dimension [U] = lehgth"^. 

Because the aim is to investigate dynamical behavior of our Universe in 
the presence of scalar fields, let us suppose that scalar fields are homoge- 
neous: (p^ = (p\t) and four-dimensional metric is spatially- flat Friedmann- 
Robertson- Walker one: g''^^ = —dt dt + o?{t)dx ® dx. 

For energy density and pressure follows: 

SttG \2 



^G.jip'ip^ = 47rG(p + P) 



(A3.3) 



U = 47cG{p - P) . 
The Friedmann equations for considered model are 



3 = 3//2 = 8t,Gp = ^GyW' + U , (A3.4) 

and 

H = -47rG{p + P) = • (A3.5) 

From these two equations, one obtains the following expression for the accel- 
eration parameter: 

^ -4 X \g,j(p'(p^ + 2[/ ) . (A3.6) 
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It can be easily seen that the equation of state (EoS) parameter uj = P j p 
and parameter q are hnearly connected: 

g = -i(l + 3a;). (A3.7) 

From the definition of the acceleration parameter, it follows that q is constant 
in the case of the power-law and De Sitter-like behavior: 



1 : a oc 



For example, q = —0.5 during the matter dominated (MD) stage where 
s = 2/3. 

Because the minisuperspace metric Gij is fiat, the scalar field equations 
are: 

+ 3Hip' + G'^^ = . (A3.9) 
For the action (IA3.ip . the corresponding Hamiltonian is 

where 



n = —G'^P.Pj + —U , (A3.10) 



i^. = ^^..^- (A3.11) 

are the canonical momenta and equations of motion have also the canonical 
form 

It can be easily seen that the latter equation (for Pi) is equivalent to the eq. 

(mm . 

Thus, the Friedmann equations together with the scalar field equations 
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can be replaced by the system of the first order ODEs: 



^ dU ^ 

a =aH, (A3.15) 

H =--H^ 
a 

4 X iG.jip'cp' + 2U] - (A3. 16) 



6 V ' 2 / 
with Eq. (IA3.4P considered in the form of the initial conditions: 



(A3.17) 



One can make these equations dimensionless: 

= rip 

Mpidt Ml^a^ 



= ^-G-P,- (A3.18) 
dt ^ ^ ^ 

dP, a^Mlid{U/Mli) 



Mpidt Stt d^' 

dP dU 



(A3.19) 



dt Stt dip^ 

That is to say the time t is measured in the Planck times tpi, the scale factor 
a is measured in the Planck lengths Lpi and the potential U is measured in 
the Mpi units. 

The system of dimensionless first order ODEs together with the initial con- 
dition (IA3.17P is used for numerical calculation of the dynamics of considered 
models with the help of a Mathematica package [ |120] . 



Appendix 4. Bouncing inflation in the Brans-Dicke frame 

According to Eq. (12.61 ). the four-dimensional FRW metrics in the Einstein 
frame (14.621) and in the Brans-Dicke frame are related as follows: 

- dt0dt + a^{t)dx 0dx = f [-dt dt + d^{t)dx dx] , (A4.1) 
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where f = X + 1 > and X is parameterized by Eq. (14.641 ). Therefore, for 
the synchronous times and scale factors in both frames one obtains, corre- 
spondingly: 

di = dt/^/m, (A4.2) 
-a{t) = am)/y/fm). (A4.3) 

which lead to the following equations: 

i= [ \ +1, (A4.4) 

Ji ^X{t) + 1 ^ ^ 

where the constant of integration is chosen in such a way that i{t = 1) = 1, 
and 



H{F) = §i = y'X («(*)) + 1 

at a 



1 dX , , ^ . 



Him) 

(A4.5) 



2{X{t{t)) + l) dt 

From the latter equation follows the relation between the Hubble parame- 
ters in both frames. In Fig. 14.331 the logarithms of the scale factor a{i) 
and the Hubble parameter H{t) ar depicted for the trajectory drawn in Fig. 
14.271 These pictures clearly demonstrate that in the Brans-Dicke frame the 
Universe has also asymptotical de Sitter stage when the scalaron field ap- 
proaches the branching point 6 = 27r. It is not difficult to verify that because 
X{t — > +oo) — > Xjnax and dX/dt{t — > +oo) — > 0, one obtains the following 
relation for the asymptotic values of the Hubble parameters in both frames: 

H = H^Xmax + 1 • (A4.6) 



Ill 



LITERATURE 

[1] Green M.B. Superstring theory / Green M.B., Schwarz J.H., Witten E. 
/ / Cambridge University Press. - 1987. 

[2] Polchinski J. String theory / Polchinski J. // Cambridge University Press. 
- 1998. 

[3] Wang L. Cosmic Concordance and Quintessence / Wang L., Caldwell 
R.R., Ostriker J. P., Steinhardt P.J. // Astrophys. J. - 2000. - V. 530. - P. 
17-35. 

[4] Arkani-Hamed N. Stabilization of submillimeter dimensions: The new 
guise of the hierarchy problem / Arkani-Hamed N., Dimopoulos S., March- 
Russell G.J. // Phys. Rev. D. - 2000. - V. 63. - 064020, 13 p. 

[5] Banks T. Constraints on theories with large extra dimensions / Banks T., 
Dine M., Nelson A.E. // JHEP. - 1999. - V. 014. - 9906, 32 p. 

[6] Arkani-Hamed N. Rapid asymmetric inflation and early cosmology in the- 
ories with sub-millimeter dimensions / Arkani-Hamed N., Dimopoulos S., 
Kaloper N., March-Russell J. // Nucl. Phys. B. - 2000. - V. 567. - P. 
189-228. 

[7] Carroll S.M. Classical stabilization of homogeneous extra dimensions / 
Carroll S.M., Geddes J., Hoffman M.B., Wald R.M. // Phys. Rev. D. - 
2002. V. 66 - 024036, 25 p. 

[8] Geddes J. The Collapse of Large Extra Dimensions / Geddes J. // Phys. 
Rev. D. - 2002. - V. 65. - 104015, 12 p. 

[9] Demir D.A. Flavordynamics with Conformal Matter and Gauge Theories 
on Compact Hyperbolic Manifolds in Extra Dimensions / Demir D.A., 
Shifman M. // Phys. Rev. D. - 2002. - Vo. 65. - 104002, 16 p. 



112 



[10] Nasri S. Radion stabilization in compact hyperbolic extra dimensions / 
Nasri S., Silva P.J., Starkman G.D., Trodden M. // Phys. Rev. D. - 2002, 
V. 66. - 045029, 16 p. 

[11] Perivolaropoulos L. Cosmological effect of radion oscillations / 
Perivolaropoulos L., Sourdis C. // Phys. Rev. D. - 2002. - V. 66. - 084018, 
7 p. 

[12] Capozziello S. A Bird's Eye View of f(R)-Gravity / Capozziello S., De 
Laurentis M., Faraoni V. // The Open Astronomy J. - 2009. - V. 3. - 
2009. - P. 49-72. 

[13] Spergel D.N. Wilkinson Microwave Anisotropy Probe (WMAP) Three 
Year Results: Imphcations for Cosmology / Spergel D.N. / / Astrophys. J. 
Suppl. - 2007. - V. 170. - 377, 91 p. 

[14] Riess A.G. Type la Supernova Discoveries at z > 1 From the Hubble 
Space Telescope: Evidence for Past Deceleration and Constraints on Dark 
Energy Evolution / Riess A.G. // Astrophys. J. - 2004. - V. 607. - P. 
665-687. 

[15] Astier P. The Supernova Legacy Survey / Astier P. // Astron. Astrophys. 
- 2006. - V. 447. - P. 31-48. 

[16] Eisenstein D.J. Sloan Digital Sky Survey (SDSS) / Eisenstein D.J. // 
Astrophys. J. - 2005. - V. 633. - P. 560-574. 

[17] Sotiriou T.P. f(R) theories of gravity / Sotiriou T.P., Faraoni V. // Rev. 
Mod. Phys. - 2010. - V. 82. - P. 451-497. 

[18] Wald R.M. General Relativity / Wald R.M. // University of Chicago 
Press. - 1984. 

[19] Guth A.H. Inflationary universe: A possible solution to the horizon and 
flatness problems / Guth A.H. // Phys. Rev. D. - 1981. - V. 23. - P. 
347-356. 



113 



[20] Kolb E.W. The Early Universe / Kolb E.W., Turner M.S. // Addison- 
Wesley. - 1992. 

[21] Linde A. Particle Physics and Inflationary Cosmology / Linde A. // 
Harwood Academic Publishers. - 1990. 

[22] Nojiri S. Unifled cosmic history in modifled gravity: from F(R) theory 
to Lorentz non-invariant models / Nojiri S., Odintsov S.D. // E Print 
arXiv.org: 1011.0544v3. 

[23] Carrofl S.M. The Cosmological Constant / Carroll S.M. // Living Rev. 
Rel. - 2001. - V. 4, 50 p. 

[24] Bondi H. Cosmology / Bondi H. // Cambridge University Press: Cam- 
bridge. - 1952. 

[25] Brans C.H. Mach's Principle and a Relativistic Theory of Gravitation / 
Brans C.H., Dicke R.H. // Phys Rev. - 1961. - V. 124. - P. 92-935. 

[26] Capozziello S. Noether Symmetries in Cosmology / Capozziello S., de 
Ritis R., Rubano C, Scudellaro P. // Riv Nuovo Cimento. - 1996. - V. 19. 
114 p. 

[27] Sciama D.W. On the Origin of Inertia / Sciama D.W. // Mon Not R 
Ast Soc. - 1953. - V. 113. - P. 34-42. 

[28] Birrell N.D. Quantum Fields in Curved Space / Birrell N.D., Davies 
P.C.W. // Cambridge University Press. - 1982. 

[29] Nojiri S. Where new gravitational physics comes from: M-theory? / 
Nojiri S., Odintsov S.D. // Phys. Lett. B. - 2003. - V. 576. - P. 5-11. 

[30] Woodard R.P. Leading Log Solution for Inflationary Yukawa / Woodard 
R.P. // Lect. Notes Phys. - 2007. - V. 720. - 044019, 35 p. 

[31] Nojiri S. Introduction to Modifled Gravity and Gravitational Alternative 
for Dark Energy / Nojiri S., Odintsov S.D. // Int. J. Geom. Meth. Mod. 
Phys. - 2007. - V. 4. - P. 115-146. 



114 



[32] Sotiriou T.P. f(R) Theories Of Gravity / Sotiriou T.P., Faraoni V. // 
Rev. Mod. Phys. - 2010. - V. 82. - P. 451-497. 

[33] Kolb E.W. Time variation of fundamental constants, primordial nucle- 
osynthesis, and the size of extra dimensions / Kolb E.W., Perry M.J., 
Walker T.P. // Phys. Rev. D. - 1986. - V. 33. - P. 869-871. 

[34] Gunther U. Gravitational excitons from extra dimensions / Gunther U., 
Zhuk k. II Phys. Rev. D. - 1997. - V. 56. - P. 6391-6402. 

[35] Gunther U. Stabilization of internal spaces in multidimensional cosmol- 
ogy / Gunther U., Zhuk A. // Phys. Rev. D. - 2000. - V. 61. - 124001, 
lip. 

[36] Gunther U. Asymptotical AdS from non-linear gravitational models with 
stabilized extra dimensions / Gunther U., Moniz P., Zhuk A. // Phys. Rev. 
D. - 2002. - V. 66. - 044014, 12 p. 

[37] Gunther U. AdS and stabilized extra dimensions in multidimensional 
gravitational models with nonlinear scalar curvature terms 1/i? and / 
Gunther U., Zhuk A., Bezerra V.B., Romero C. // Class. Quant. Grav. - 
2005. - V. 22. - 3135, 28 p. 

[38] Kerner R. Cosmology without singularity and non-linear gravitational 
lagrangians / Kerner R. // GRG. - 1982. - V. 14. - P. 453-469. 

[39] Barrow J.D. The stability of general relativistic cosmological theory / 
Barrow J.D., Ottewill A.C. // J. Phys. A. - 1983. - V. 10. - P. 2757-2776. 

[40] Duruisseu J. P. Non-Einsteinian lagrangians assuring non-singular cos- 
mological solutions / Duruisseu J. P., Kerner R. // GRG. - 1983. - V. 15. 
- P. 797-807. 

[41] Whitt B. Fourh-order gravity as general relativity plus matter / Whitt 
B. // Phys. Lett. B. - 1984. - V. 145. - P. 176-178. 



115 



[42] Barrow J.D. Inflation and the conformal structure of higher-order gravity 
theories / Barrow J.D., Cotsakis S. // Phys. Lett. B. - 1988. - V. 214. - 
P. 515-518. 

[43] Maeda K. Chaotic inflationary scenario of the Universe with a nonmini- 
mally coupled "inflaton" fleld / Maeda K., Stein-Schabes J. A., Futamase 
T. // Phys. Rev. D. - 1989. - V. 39. - P. 2848-2853. 

[44] Magnano G. On Physical Equivalence between Nonlinear Gravity The- 
ories / Magnano G., Sokolowski L.M. // Phys. Rev. D. - 1994. - V. 50. - 
P. 5039-5059. 

[45] Wands D. Extended Gravity Theories and the Einstein-Hilbert Action / 
Wands D. // CQG. - 1994. - V. 11. - P. 269-279. 

[46] Elhs J. Topological i?^ Inflation / Elhs J., Kaloper N., Olive K.A., 
Yokoyama J. // Phys. Rev. D. - 1999. - V. 59. - 103503, 54 p. 

[47] Nojiri S. Dark energy, inflation and dark matter from modifled F(R) 
gravity / Nojiri S., Odintsov S.D. // E Print arXiv.org: 0807.0685. 

[48] Witt B. Fourth-order gravity as general relativity plus matter / Witt B. 
// Phys. Lett. B. - 1984. - V. 145. - P. 176-178. 

[49] Muller V. The stability of the de Sitter space-time in fourth order gravity 
/ Muller v., Schmidt H.J., Starobinsky A.A. // Phys. Lett. B. - 1988. - 
V. 202. - P. 198-200. 

[50] Maeda K. Towards the Einstein-Hilbert action via conformal transfor- 
mation // Phys. Rev. D. - 1989. - V. 39. - P. 3159-3162. 

[51] Kaneda S. Fourth-order gravity as the inflationary model revisited / 
Kaneda S., Ketov S.V., Watanabe N. // Mod. Phys. Lett. A. - 2010. - V. 
25. - P. 2753-2762. 

[52] Kaneda S. Slow-roll inflation in (i? + R^) gravity / Kaneda S., Ketov 
S.V., Watanabe N. // Class. Quantum Grav. - 2010. - V. 27. - 145016, 12 
P- 



116 



[53] Starobinsky A. A. A new type of isotropic cosmological model without 
singularity / Starobinsky A. A. // Phys. Lett. B. - 1980. - V. 91. - 99, 4 p. 

[54] Saidov T. 1/R multidimensional gravity with form-fields: Stabilization 
of extra dimensions, cosmic acceleration, and domain walls / Saidov T., 
Zhuk. k. II Phys. Rev. D. - 2007. - V. 75. - 084037, 10 p. 

[55] Carroll S.M. The Cosmology of Generalized Modified Gravity Models / 
Carroll S.M., De Felice A., Duvvuri V., Easson D.A., Trodden M., Turner 
M.S. // Phys.Rev. D. - 2005. - V. 71. - 063513, 27 p. 

[56] Elhs J. Topological i?^ Infiation / Elhs J., Kaloper N., Olive K.A., 
Yokoyama J. // J. Phys. Rev. D. - 1999. - V. 59. - 103503, 25 p. 

[57] Gunther U. Asymptotical AdS from non-linear gravitational models with 
stabilized extra dimensions / Gunther U., Moniz P., Zhuk A. // Phys. Rev. 
D. - 2002. - V. 66. - 044014, 12 p. 

[58] Gunther U. Nonlinear multidimensional cosmological models with form 
fields: stabilization of extra dimensions and the cosmological constant prob- 
lem / Gunther U., Moniz P., Zhuk A. // Phys. Rev. D. - 2003. - V. 68. - 
044010, 21 p. 

[59] Saidov T. AdS Nonlinear Curvature- Squared and Curvature-Quartic 
Multidimensional (D=8) Gravitational Models with Stabilized Extra Di- 
mensions / Saidov T., Zhuk. A. // Gravitation and Cosmology. - 2006. - 
V. 12. - P. 253-261. 

[60] Bronnikov K.A. Abilities of multidimensional gravity / Bronnikov K.A., 
Rubin S.G. // Grav. and Cosmol. - 2007. - V. 13. - 253, 6 p. 

[61] Saidov T. Problem of infiation in nonlinear multidimensional cosmologi- 
cal models / Saidov T., Zhuk A. // Phys. Rev. D. - 2009. - V. 79. - 024025, 
18 p. 

[62] Arkani-Hamed N. Phenomenology, Astrophysics and Cosmology of The- 
ories with Sub-Millimeter Dimensions and TeV Scale Quantum Gravity / 



117 



Arkani-Hamed N., Dimopoulos S., Dvali G. // Phys. Rev. D. - 1999. - V. 
59. - 086004, 52 p. 

[63] Gunther U. Multidimensional cosmological models: cosmological and 
astrophysical implications and constraints / Gunther U., Starobinsky A., 
Zhuk A. // Phys. Rev. D. - 2004. - V. 69. - 044003, 16 p. 

[64] Srianand R. Aracil B. Limits on the time variation of the electromagnetic 
fine-structure constant in the low energy limit from absorption lines in the 
spectra of distant quasars / Srianand R., Chand H., Petitjean P., Aracil B. 
// Phys. Rev. Lett. - 2004. - V. 92. - 121302, 4 p. 

[65] Hannestad S. Possible Constraints on the Time Variation of the Fine 
Structure Constant from Cosmic Microwave Background Data / Hannestad 
S. // Phys. Rev. D. - 1999. - V. 60. - 023515, 5 p. 

[66] Kaplinghat M. Constraining Variations in the Fine-Structure Constant 
with the Cosmic Microwave Background / Kaplinghat M., Scherrer R.J., 
Turner M.S. // Phys. Rev. D. - 1999. - V. 60. - 023516, 9 p. 

[67] Webb J.K. Further Evidence for Cosmological Evolution of the Fine 
Structure Constant / Webb J.K. // Phys. Rev. Lett. - 2001. - V. 87. - 
091301, 5 p. 

[68] Ivashchuk V.D. On Wheeler-De Witt equation in multidimensional cos- 
mology / Ivashchuk V.D., Melnikov V.N., Zhuk A.I. / / Nuovo Cimento B. 
- 1989. - V. 104. - P. 575-582. 

[69] Rainer M. Tensor-multi-scalar theories from multidimensional cosmology 
/ Rainer M., Zhuk A. // Phys. Rev. D. - 1996. - V. 54. - P. 6186-6192. 

[70] Capozziello S. Quintessence without scalar field / Capozziello S., Carloni 
S., Troisi A. // Recent Res. Dev. Astron. Astrophys. - 2003. - V. 1. - 625, 
46 p. 



118 



[71] Carrol S.M. Is Cosmic Speed-Up Due to New Gravitational Physics? / 
arrol S.M., Duvvuri V., Trodden M., Turner M.S. // Phys. Rev. D. - 2004. 

- V. 70. - 043528, 4 p. 

[72] Vollick D.N. 1/R Curvature Corrections as the Source of the Cosmologi- 
cal Acceleration / Volhck D.N. // Phys. Rev. D. - 2003. - V. 68. - 063510, 
7 p. 

[73] Dick R. On the Newtonian Limit in Gravity Models with Inverse Powers 
of R / Dick R. // Gen. Rel. Grav. - 2004. - V. 36. - 217, 8 p. 

[74] Dolgov A.D. Can modified gravity explain accelerated cosmic expansion? 
/ Dolgov A.D., Kawasaki M. // Phys. Lett. B. - 2003. - V. 573. - P. 1-4. 

[75] Nojiri S. Modified gravity with negative and positive powers of the cur- 
vature: unification of the infiation and of the cosmic acceleration / Nojiri 
S., Odintsov S.D. // Phys. Rev. D. - 2003. - V. 68. - 123512, 23 p. 

[76] Nojiri S. Modified gravity with InR terms and cosmic acceleration / 
Nojiri S., Odintsov S.D.// Gen. Rel. Grav. - 2004. - V. 36. - P. 1765- 
1780. 

[77] Chiba T. 1/R gravity and Scalar-Tensor Gravity / Chiba T.// Phys. 
Lett. B. - 2003. - V. 575, 4 p. 

[78] Meng X. Modified Friedmann Equations in R~^-Modified Gravity / Meng 
X., Wang P.// Class. Quant. Grav. - 2003. - V. 20. - P. 4949-4962. 

[79] Gunther U. Multidimensional Cosmology and Asymptotical AdS / Gun- 
ther U., Moniz P., Zhuk A. // Astrophysics and Space Science. - 2003. - 
V. 283. - P. 679-684. 

[80] Freund P.G.O. Dynamics of dimensional reduction / Freund P.G.O., Ru- 
bin M.A. // Phys. Lett. B. - 1980. - V. 97. - P. 233-235. 

[81] Gunther U. Multidimensional perfect fiuid cosmology with stable com- 
pactified internal dimensions / Gunther U., Zhuk A. // Class. Quant. Grav. 

- 1998. - V. 15. - 2025, 10 p. 



119 



[82] Ivashchuk V.D. Billiard Representation for Multidimensional Cosmology 
with Multicomponent Perfect Fluid near the Singularity / Ivashchuk V.D., 
Melnikov V.N. // Class. Quant. Grav. - 1995. - V. 12. - 809, 16 p. 

[83] Zhuk A. Integrable scalar field multi-dimensional cosmologies / Zhuk A. 
// Class. Quant. Grav. - 1996. - V. 13. - 2163, 17 p. 

[84] Bousso R. Quantization of Four-form Fluxes and Dynamical Neutraliza- 
tion of the Cosmological Constant / Bousso R., Polchinski J. // JHEP. - 
2000. - V. 006. - 0006, 30 p. 

[85] Schwartz-Perlov D. Probabilities in the Bousso-Polchinski multiverse / 
Schwartz-Perlov D., Vilenkin A. // JCAP. - 2006. - V. 010. - 0606, 22 p. 

[86] Zhuk A. Conventional cosmology from multidimensional models / Zhuk 
A. // E Print arXiv.org: hep-th/0609126. 

[87] Linde A. Monopoles as big as a universe / Linde A. // Phys. Lett. B. - 
1994. - V. 327. - P. 208-213. 

[88] Vilenkin A. Topological inflation / Vilenkin A. // Phys. Rev. Lett. - 
1994. - V. 72. - P. 3137-3140. 

[89] Sakai N. Dynamics of Topological Defects and Inflation / Sakai N., 
Shinkai H., Tachizawa T., Maeda K. // Phys.Rev. D. - 1996. - V. 53. 
- 655, 15 p. 

[90] Hindawi A. Higher-Derivative Gravitation and a New Mechanism for 
Supersymmetry Breaking in Four-Dimensions / Hindawi A., Ovrut B.A., 
Waldram D. // Prog. Theor. Phys. Suppl. - 1996. - V. 123. - P. 397-410. 

[91] Kiritsis E. On threshold corrections in IIB string theory and (p,q) 
string instantons / Kiritsis E., Pioline B. // Nucl. Phys. B. - 1997. - V. 
508. - P. 509-534. 

[92] Green M.B. One loop in eleven dimensions / Green M.B., Gutperle M., 
Vanhove P. // Phys. Lett. B. - 1997. - V. 409. - P. 177-184. 



120 



[93] Bachas CP. Curvature terms in D-brane actions and their M-theory 
origin / Bachas CP., Bain P., Green M.B. // JHEP. - 1999. - V. Oil. - 
9905, 32 p. 

[94] Nojiri S. Finite gravitational action for higher derivative and stringy 
gravities / Nojiri S., Odintsov S. // Phys. Rev. D. - 2000. - V. 62. - 
064018, 21 p. 

[95] Zhuk A. Restrictions on dilatonic brane- world models / Zhuk A. // Int. 
Journ. Mod. Phys. D. - 2002. - V. 11. - 1399, 10 p. 

[96] Baukh V. Sp-brane accelerating cosmologies / Baukh V., Zhuk A. // 
Phys. Rev. D. - 2006. - V. 73. - 104016, 17 p. 

[97] Rubakov V.A. Large and infinite extra dimensions / Rubakov V.A. // 
Phys. Usp. - 2001. - V. 44. - P. 871-893. 

[98] Appelquist T. Bounds on Universal Extra Dimensions / Appelquist T. , 
Cheng H.C, Dobrescu B.A. // Phys.Rev. D. - 2001. - V. 64. - 035002, 22 

P- 

[99] Ponton E. Casimir Energy and Radius Stabilization in Five and Six 
Dimensional Orbifolds / Ponton E., Poppitz E. // JHEP. - 2011. - V. 019. 
- 0106, 38 p. 

[100] Cheng H.C. Schmaltz M. Radiative Corrections to Kaluza-Klein Masses 
/ Cheng H.C, Matchev K.T., // Phys. Rev. D. - 2002. - V. 66. - 036005, 
33 p. 

[101] Servant G. Is the Lightest Kaluza-Klein Particle a Viable Dark Matter 
Candidate? / Servant C, Tait T.M.P. // Nucl.Phys. B. - 2003. - V. 650. 
-391, 32 p. 

[102] Bergstrom L. Gamma Rays from Kaluza-Klein Dark Matter / 
Bergstrom L., Bringmann T., Eriksson M., Gustafsson M. // Phys. Rev. 
Lett. - 2005. - V. 94. - 131301, 14 p. 



121 



[103] Townsend P.K. Four lectures on M-theory / Townsend P.K. // E Print 
arXiv.org: hep-th/9612121. 

[104] Moffat J.W. M-theory on a supermanifold / Moffat J.W. // E Print 
arXiv.org: hep-th/0111225. 

[105] Kundu A. Universal Extra Dimension / Kundu A. // E Print arXiv.org: 
0806.3815. 

[106] Ratra B. Cosmological consequences of a rolling homogeneous scalar 
field / Ratra B., Peebles P.J. // Phys. Rev. D. - 1988. - V. 37. - P. 3406- 
3427. 

[107] Elhs G.F.R. Exact scalar field cosmologies / Ellis G.F.R., Madsen M.S. 
// Class. Quant. Grav. - 1991. - V. 8. - 667, 11 p. 

[108] Ratra B. Inflation in an exponential-potential scalar field model / Ratra 
B. // Phys. Rev. D. - 1992. - V. 45. - P. 1913-1952. 

[109] Stornaiolo C. Cosmological fiuids with general equations of state / Stor- 
naiolo C. // Phys. Let. A. - 1994. - V. 189. - P. 351-360. 

[110] Gunther U. Gravitational excitons from extra dimensions / Gunther 
U., Zhuk k. II Phys. Rev. D. - 1997. - V. 56. - 6391, 20 p. 

[Ill] Arkani-Hamed N. Stabilization of Sub-Millimeter Dimensions: The 
New Guise of the Hierarchy Problem / Arkani-Hamed N., Dimopoulos S., 
March-Russeh J. // Phys.Rev. D. - 2001. - V. 63. - 064020, 30 p. 

[112] Cho Y.M. Dilatonic dark matter and unified cosmology: a new 
paradigm / Cho Y.M., Keum Y.Y. // Class. Quant. Grav. - 1998. - V. 15. 

- 907, 16 p. 

[113] Cho Y.M. Dilaton as a Dark Matter Candidate and its Detection / Cho 
Y.M., Kim J.H. // Phys. Rev. D. - 2009. - V. 79. - 023504, 23 p. 

[114] Martin J. The precision of slow-roll predictions for the CMBR 
anisotropics / Martin J., Schwarz D. // Phys. Rev. D. - 2000. - V. 62. 

- 103520, 15 p. 



122 



[115] Gunther U. Stabilization of internal spaces in multidimensional cosmol- 
ogy / Gunther U., Zhuk A. // Phys. Rev. D. - 2000. - V. 61. - 124001, 11 
P- 

[116] DeruelleN. "Detuned" f(R) gravity and dark energy / DeruelleN., Sasaki 
M., Sendouda Y. // Phys. Rev. D. - 2008. - V. 77. - 124024, 10 p. 

[117] Nakamura T.T. The spectrum of cosmological perturbations produced 
by a multi-component inflaton to second order in the slow-roll approxima- 
tion / Nakamura T.T., Stewart E.D. // Phys. Lett. B. - 1996. - V. 381. - 
P. 413-419. 

[118] Gong J.O. The Power Spectrum for a Multi-Component Inflaton to 
Second-Order Corrections in the Slow-Roll Expansion / Gong J.O., Stewart 
E.D. // Phys. Lett. B. - 2002. - V. 538. - P. 213-222. 

[119] Nibbelink S.G. Density perturbations arising from multiple field slow- 
roll infiation / Nibbelink S.G., Tent B.J.W. // E Print arXiv.org: hep- 
ph/0011325v2. 

[120] Kallosh R. SuperCosmology / Kallosh R., Prokushkin S. // E Print 
arXiv.org: hep-th/0403060. 

[121] Lyth D.H. Particle Physics Models of Inflation and the Cosmological 
Density Perturbation / Lyth D.H., Riotto A. // Phys. Rept. - 1999. - V. 
314, 155 p. 

[122] Efstathiou G. The Lyth Bound Revisited / Efstathiou G., Mack K.J. 
// JCAP. - 2005. - V. 008. - 0505, 8 p. 

[123] Peiris H.V. Primordial Black Holes, Eternal Inflation, and the Inflation- 
ary Parameter Space after WMAP5 / Peiris H.V., Easther R. // JCAP. - 
2008. - V. 024. - 0807, 24 p. 

[124] Bronnikov K.A. Abilities of multidimensional gravity / Bronnikov K.A., 
Rubin S.G. // Grav and Cosmol. - 2007. - V. 13. - 253, 6 p. 



123 



[125] Linde A. Monopoles as Big as a Universe / Linde A. // Phys. Lett. B. 

- 1994. - V. 327. - P. 208-213. 

[126] Vilenkin A. Topological Inflation / Vilenkin A. // Phys. Rev. Lett. - 
1994. - V. 72. - P. 3137-3140. 

[127] Liddle A.R. The Cold Dark Matter Density Perturbation / Liddle A.R., 
Lyth D.H. // Phys. Rept. - 1993. - V. 231, 105 p. 

[128] Barenboim G. About a (standard model) universe dominated by the 
right matter / Barenboim G., Vives O. // Phys. Rev. D. - 2009. - V. 79. 

- 033007, 15 p. 

[129] Frolov A.V. A Singularity Problem with f(R) Dark Energy / Frolov 
A.V. // Phys. Rev. Lett. - 2008. - V. 101. - 061103, 5 p. 

[130] Goleman S. Gravitational effects on and of vacuum decay / Goleman 
S., De Luccia F. // Phys. Rev. D. - 1980. - V. 21. - P. 3305-3315. 

[131] Gai Y. A Model Of Inflationary Cosmology Without Singularity / Cai 
Y., Qiu T., Xia J., Zhang X. // Phys. Rev. D. - 2009. - V. 79. - 021303, 
4 p. 

[132] Paul B.C. Singular instantons in higher derivative theories / Paul B.C., 
Mukherjee S., Tavakol R.K. // Phys. Rev. D. - 2002. - V. 65. - 064020, 
14 p. 

[133] Biswas T. Bouncing Universe in string-inspired gravity / Biswas T., 
Mazumdar A., Siegel W. // JCAP. - 2006. - V. 009. - 0603, 28 p. 

[134] Lyth D.H. Large-scale energy-density perturbations and inflation / 
Lyth D.H. // Phys. Rev. D. - 1985. - V. 31. - P. 1792-1798. 



124 



ACKNOWLEDGMENT 



First of all, I would like to thank my parents, especially my mother, who 
instilled in me a love for science. 

With my greatest pleasure I would like to thank my supervisor Prof. 
Alexander Ivanovich Zhuk, who has opened for me the multidimensional cos- 
mology, for his strong support in the research, fruitful discussions and many 
of very interesting walk tours we had around Sevastopol and Odessa. His 
capability for elegant solutions in very complex theoretical tasks will always 
amaze me. 

I would like to address my warmest gratitude to professors of the De- 
partment of Theoretical Physics ONU: Prof. V.M. Adamyan, Prof. N.P. 
Malomuzh and Prof. A.V. Zatovskiy , Dr. V.P. Olyiynik, Dr. M.Ya. Sushko 
and Dr. V.L. Kulinskiy; their lectures allowed me to enter the world of 
theoretical physics. 

Then I like to thank my Teacher of physics (deserving capital T) Vadim 
Leonidovich Manakin, whose impact on my life is impossible to overestimate. 
There are people gifted with ability to kindle a passion to knowledge and 
science. I was lucky to meet such a person - Vadim Leonidovich. 

And last, but definitely not least, I would like to acknowledge my close 
friend Volodymyr Volodymyrovych Nemertsalov for his great help in prepa- 
ration of promotion ceremony. 

Thank you all! 



125 



CURRICULUM VITAE 

Tamerlan Adamovich Saidov was born in Grozny (Chechen Repubhc) in 
October 1, 1982. In 1999 he graduated with the certificate of merit the Riche- 
heu Lyceum (Department of Chemistry and Biology) in Odessa, Ukraine. 

In 2003 he received the Bachelor degree in Theoretical Physics with Ex- 
cellence at the Faculty of Physics (Department of Theoretical Physics) of 
the Odessa I.I. Mechnikov National University (Odessa, Ukraine). The study 
of the bachelor report was focused on Navier-Stokes equations in cylindrical 
coordinates. 

In period of 2003 - 2004, Tamerlan Saidov has obtained the Master degree 
in Theoretical Physics with Excellence at the same department. The work 
of the master program was dedicated to the properties of fermions fields in 
models of charged leptons, and was performed under supervision of Dr. V.P. 
Olyiynik. 

Simultaneously with studying of physics, in 2005 Tamerlan Saidov gradu- 
ated the Faculty of Economy and Law (Department of Administrative Law) 
of Odessa I.I. Mechnikov National University (Odessa, Ukraine) and received 
Master degree in Law. The master thesis was written on the administra- 
tive right and obligations with the focus on the specifics of their origin and 
prosecuting. 

In the period of 2004 - 2008, Tamerlan Saidov was appointed as phd- 
student at the Department of Theoretical Physics of the Odessa I.I. Mech- 
nikov National University (Odessa, Ukraine), where under the supervision 
of Prof. A.I. Zhuk, he has undertaken the theoretical study on multidimen- 
sional cosmology: the compactification of extra dimensions and accelerated 
expansion of the Universe. The results of this work are described in this 
Thesis. 

Starting from 2007 to 2011, he did his second phd-research, in the Depart- 
ment of Applied Physics of Eindhoven University of Technology (Eindhoven, 
The Netherlands). Where under the supervision of Dr. L. Pel, the crystal- 
lization and ion transport of salts was studied in porous media by Nuclear 



126 



Magnetic Resonance technique. 

From 2011, Tamerlan Saidov is working on the cancer detection by con- 
trast ultrasound imaging in the Signals Processing Systems group at the 
Department of Electrical Engineering of Eindhoven University of Technology 
(Eindhoven, The Netherlands). 



